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ABSTRACT 

The  paper  describes  the  development  of  an  aeroelastic  model  for  the  prediction  of  empennage 
dynamic  response  due  to  vortex  induced  buffet  loads.  The  multidisciplinary  problem  of  tail 
buffeting  is  solved  accurately  in  the  time  domain  using  an  unsteady  vortex  model  for 
prediction  of  aerodynamic  loads  and  coupled  aeroelastic  equations  for  the  bending  and 
torsional  deflections  of  the  tail  which  are  resolved  using  the  Galerkin  method.  A  dynamic 
aeroelastic  analysis  of  empennage  buffet  is  performed  for  a  generic  delta  wing,  twin 
vertical-tail  configuration  at  high  angles  of  attack.  Results  of  computations  are  validated  using 
experimental  and  flight  test  data.  The  present  research  provides  advances  in  the  predictive 
capability  and  our  understanding  of  aircraft  empennage  buffet. 
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Vertical  Tail  Dynamic  Response  in  Vortex 
Breakdown  Flow 


Executive  Summary 

The  airframes  of  high  performance  aircraft,  such  as  the  F/A-18,  have  suffered  from  an 
aeroelastic  tail  buffet  problem  for  many  years.  This  problem  is  inherent  to  vortical 
flows  used  to  generate  lift  at  high  angles  of  attack  as  they  tend  to  break  down  causing 
severe  empennage  dynamic  loading  and  premature  fatigue  failures. 

Better  understanding  of  the  empennage  buffeting  problem  is  required  for  development 
of  reliable  fatigue  usage  monitoring  systems  and  for  the  fleet  management  of  aircraft. 
The  challenges  associated  with  computational  simulation  of  empennage  buffet  vary 
from  prediction  of  the  nonlinear  separated  vortical  flows  about  complex  configurations 
to  the  coupled  interaction  between  the  flow  and  the  dynamic  response  of  the  tail 
structure. 

The  work  describes  the  development  and  validation  of  an  aeroelastic  model  for  the 
prediction  of  the  F /  A-18  empennage  buffet  due  to  bursting  LEX  vortices  under  various 
flight  conditions.  The  multidisciplinary  problem  of  tail  buffeting  is  solved  accurately  in 
the  time  domain  using  an  unsteady  vortex  model  for  prediction  of  aerodynamic  loads 
and  coupled  aeroelastic  equations  for  the  bending  and  torsional  deflections  of  the  tail 
which  are  resolved  using  the  Galerkin  method. 

A  dynamic  aeroelastic  analysis  of  empennage  buffet  is  performed  for  a  generic  delta 
wing,  twin  vertical-tail  configuration  at  high  angles  of  attack.  Results  of  flow 
simulation  indicated  that  the  aeroelastic  model  is  able  to  predict  major  unsteady 
features  of  the  vortex  induced  buffet  loads  and  the  resulting  coupled  fluid-structure 
interaction. 

The  present  research  provides  advances  in  the  predictive  capability  and  our 
understanding  of  aircraft  empennage  buffet.  Results  of  the  work  contribute  to  DSTO's 
existing  body  of  knowledge  on  vortex  breakdown  and  can  be  used  in  assessing  the 
buffet  environments  of  current  and  future  generations  of  fighter  aircraft. 
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Nomenclature 


[/J  Identity  matrix 

\k\  Stiffness  matrix 

[m]  Mass  matrix 

A  Aerodynamic  loading  vector 

Aj  Panel  area 

Ad  Dynamic  component  of  aerodynamic  loading  vector 
Af  Tail  planform  area 

As  Static  component  of  aerodynamic  loading  vector 

Cp  Tail  root  chord 

CM  Bending  moment  coefficient 

C'M  Root-mean-square  value  of  bending  moment  coefficient 
CN  Normal  force  coefficient 

C'N  Root-mean-square  value  of  normal  force  coefficient 
Cp  Pressure  coefficient 

C'p  Root-mean-square  value  of  pressure  coefficient 

El  Distributed  bending  stiffness 

€j  Local  truncation  error 

Ek  Kinetic  energy  of  the  system 

EP  Potential  energy  of  the  system 

F  Static  loading  vector  due  to  gravity 

f  State  variable  function 

f,  Frmction  of  surface  geometry  and  flow  velocity 

fp  Pressure  function 

f„  Cut-off  function 

GJ  Distributed  torsional  stiffness 

g  Acceleration  due  to  gravity 

I  Number  of  bending  modes 

Iq  Mass  moment  of  inertia  about  the  elastic  axis 

Icg  Mass  moment  of  inertia  about  the  inertial  axis 

J  Total  number  of  modes 

L  Vertical  tail  length 

lc  Characteristic  length 

/,  Distance  from  the  tail  root  to  the  centre  of  the  panel 

M  Distributed  aerodynamic  pitch  moment 

m  Distributed  mass 

Ma  Dynamic  component  of  distributed  aerodynamic  pitch  moment 
Ms  Static  component  of  distr  ibuted  aerodynamic  pitch  moment 

N  Distributed  aerodynamic  normal  force 

n  Unit  normal  vector 

Ad  Dynamic  component  of  distributed  aerodynamic  force 
Ns  Static  component  of  distributed  aerodynamic  force 

p  Flow  pressure 

P'  Root-mean-square  value  of  pressure 

Px  Flow  pressure  at  infinity 

Q  Velocity  singular  kernel 

qi  Generalised  coordinate  for  bending  mode 


(]j  Generalised  coordinate  for  torsion  mode 

qx  Free-stream  dynamic  pressure 

r  Flow  field  coordinate  vector 

rp  Separation  line  coordinate  vector 

rs  Stagnation  point  coordinate  vector 

r  Vortex  wake  coordinate  vector 

°P 

r0  Body  surface  coordinate  vector 

s  Parameter  along  the  curve 

t  Time 

u  Flow  velocity 

U  Body  linear  velocity 

Uo  Flow  field  characteristic  velocity 

V  Induced  flow  velocity 

Voc  Flow  velocity  at  infinity 

w  Dimensionless  flow  velocity 

W  Relative  flow  velocity 

W A  Work  performed  by  the  aerodynamic  loading 

Xcc,  Coordinate  of  inertial  axis 

xc  Coordinate  of  elastic  axis 

Xo  Local  distance  between  the  elastic  and  inertial  axis 

Y  Bending  deflection 

Ya  Dynamic  component  of  bending  deflection 

7S  Static  component  of  bending  deflection 

z  Distance  from  tire  fixed  support  along  the  tail  elastic  axis 

(3  Angle  between  tire  vertical  tail  and  delta  wing  plane 

v  Flow  viscosity 

V| /,  Bending  mode  shape 

V|/j  Torsion  mode  shape 

0  Torsional  deflection 

9d  Dynamic  component  of  torsional  deflection 

9S  Static  component  of  torsional  deflection 

<)>  Dimensionless  velocity  potential 

T  Circulation  of  a  vortex  filament 

y  Vorticity  sheet  strength 

d  Vortex  core  size 

//  Time  derivative  of  generalised  coordinate 

k  Point  on  a  body  surface 

ju  Variable  point  on  a  body  surface 

C  Smoothing  function 

a  Distributed  vorticity 

a  Separated  vorticity  sheet 

p  Flow  density 

O  Velocity  potential 

fl  Body  angular  velocity 

co  Flow  vorticity 
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1.  Introduction 


For  modem  combat  aircraft,  the  ability  to  fly  and  manoeuvre  at  high  angles  of  attack 
and  at  high  loading  conditions  gives  tactical  advantages.  For  the  F/A-18  aircraft  the 
manoeuvrability  at  very  high  angles  of  attack  is  achieved  through  a  combination  of  the 
wing  root  leading  edge  extensions  (LEXs)  and  the  placement  of  twin  vertical  tails. 

A  pair  of  large  vortices  generated  by  highly  swept  wing  root  LEXs  contribute  to 
enhanced  vortex  lift  by  developing  high  suction  areas  over  the  wings.  The  F/A-18 
aircraft  is  designed  to  utilise  unique  characteristics  of  these  vortex  structures  which  are 
prime  contributors  to  the  aerodynamics  of  an  aircraft  during  high  angle  of  attack  flight. 
The  twin  vertical  tails  of  the  F/A-18  aircraft  are  placed  in  closer  proximity  to  the  LEX 
vortex  flowfield  to  take  full  advantage  of  the  concentrated  energy  contained  in  the  LEX 
vortices  to  provide  the  directional  stability  and  control  necessary  for  high  angle  of 
attack  manoeuvrability. 

However,  the  characteristics  of  these  LEX  vortices  are  highly  sensitive  to  flight 
conditions  and  there  is  a  price  for  the  enhancement  in  lift.  At  higher  angles  of  attack, 
the  LEX  vortices  are  unable  to  withstand  large  adverse  pressure  gradients  present  in 
the  flowfield  over  the  wings  and  the  initially  stable  vortices  tend  to  break  down  prior 
to  reaching  the  vertical  tails.  Flow  reversal  in  the  vortex  cores  occurs  resulting  in  a 
highly  turbulent  flow  that  impinges  on  the  tail  surfaces  causing  their  dynamic  loading 
and  vibrations  of  substantial  magnitudes.  Here,  the  appearance  of  vortex-induced 
unsteady  loads  exciting  the  tail  structure  is  commonly  known  as  buffet  while  forced 
response  of  the  tail  surfaces  to  dynamic  loading  is  often  referred  to  as  buffeting. 

While  LEX  vortices  provide  enhanced  lift  and  maneuverability,  their  breakdown  can 
cause  substantial  fatigue  damage  to  the  vertical  tails  of  the  F/A-18  aircraft  due  to 
unsteady  buffet  loading  when  the  dominant  frequency  of  vortex  breakdown  flow  is 
close  to  the  resonant  structural  frequencies  of  the  tail  structure.  Based  on  original 
designs  and  service  usage,  the  aft  fuselage  and  empennage  of  the  F/A-18  aircraft 
experienced  cracks  caused  by  vortex-induced  tail  buffet  at  a  relatively  low  number  of 
flight  hours,  and  significant  modifications  of  the  empennage  structure  were  required. 

To  alleviate  the  buffet  loads  on  the  F/A-18  vertical  tails,  the  LEX  fence  was  developed 
by  McDonnell  Douglas  Corporation  using  extensive  wind  tunnel  and  flight  tests. 
Although  a  significant  reduction  of  the  F/A-18  empennage  buffet  loads  has  been 
obtained  as  a  result  of  using  LEX  fences,  the  search  for  more  effective  solutions  is  still 
underway  for  future  generation  fighter  aircraft  as  their  design  is  still  focusing  on  high 
angle  of  attack  manoeuvrability  at  high  loading  conditions,  renewing  interest  in  the 
empennage  buffeting  problem.  Our  ability  to  understand  and  predict  empennage 
buffeting  may  not  only  help  to  alleviate  inherent  problems  associated  with  this  flow 
phenomenon  but  may  also  improve  aircraft  controllability  and  maneuverability. 
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The  technical  challenges  associated  with  die  solution  of  the  empennage  buffet  problem 
vary  from  prediction  of  the  nonlinear  separated  vortical  flows  about  complex 
configurations  to  the  coupled  interaction  between  the  flow  and  the  response  motion  of 
die  tail  structure.  Aldiough  aeroelastic  models  based  on  Euler  and  Navier-Stokes  codes 
can  be  applied  for  the  simulation  of  F/A-18  vertical  tail  buffeting,  they  require 
enormous  computing  resources  to  perform  parametric  studies.  Even  the  single 
computation  of  a  time-accurate  solution  of  F/A-18  tail  buffet  at  certain  flight 
conditions  using  die  above  codes  is  at  the  very  limit  of  the  capabilities  of  modem 
supercomputers.  An  alternative  to  computationally  intensive  Euler  and  Navier-Stokes 
codes  can  be  vortex-based  methods  which  provide  economy  in  computations  by 
concentrating  their  efforts  in  die  areas  of  high  vorticity  gradients.  The  vortex  methods 
are  particularly  suitable  for  simulation  of  vorticity  dominated  regions,  such  as  flow 
past  a  LEX. 

The  purpose  of  this  research  is  to  provide  an  improvement  in  aeroelastic  analysis 
capabilities  tiirough  die  application  of  the  unsteady  vortex  model  to  the  aeroelastic 
formulation.  This  includes  the  modeling  of  leading  edge  vortex  development  and 
breakdown,  and  the  interaction  of  a  burst  vortex  flowfield  with  flexible  vertical  tails. 
Furtiier  improvements  such  as  structural  nonlinearities  can  be  easily  incorporated  in 
die  aeroelastic  scheme.  However,  the  emphasis  in  this  research  is  placed  on  the 
aerodynamic  model  and  die  solution  strategy. 


2.  Aeroelastic  Phenomenon  of  Tail  Buffeting 

The  airframes  of  high  performance  aircraft,  such  as  the  F/A-18,  have  suffered  from  the 
aeroelastic  tail-buffeting  problem  due  to  vortex  breakdown  for  many  years.  This 
problem  is  inherent  in  any  aircraft  design  which  relies  on  the  generation  of  vortex  lift 
for  high  angle  of  attack  capabilities.  Considerable  advances  have  been  made  during  the 
past  decade  in  the  predictive  capability  and  basic  understanding  of  tail  buffet. 
However,  tiiis  aeroelastic  phenomenon  still  remains  one  of  the  most  challenging 
fundamental  research  problems  and  die  search  for  more  effective  solutions  to  the 
buffeting  problem  is  still  under  way  for  current  and  future  generations  of  fighter 
aircraft. 

Tire  term  'tail  buffeting'  was  first  introduced  after  an  accident  resulting  in  the  loss  of  a 
small  transport  airplane  at  Meopham,  England  in  1930  [1].  After  detailed  investigation 
it  was  concluded  tiiat  the  most  probable  cause  of  the  accident  was  the  horizontal  tail 
structural  failure  due  to  its  intense  forced  vibrations  or  'empennage  buffet'  caused  by 
turbulent  flow  separated  from  wings  and  fuselage.  It  was  suggested  that  the  aircraft, 
while  in  level  tiight,  entered  the  region  of  a  strong  upward  gust  that  resulted  in  a  sharp 
increase  of  angle  of  attack  and  subsequent  dow  separation  over  die  wing  causing  tail 
buffet. 
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Due  to  the  nature  of  its  geometry,  tail  buffeting  caused  by  separating  flow  over  the 
wing  was  the  major  concern  for  an  early  fighter-type  aircraft  where  significant  buffet 
usually  occurred  for  moderate  angles  of  attack  at  high  subsonic  speeds.  The  critical 
combination  of  buffeting  and  maneuver  loads  often  occurred  while  attaining  high 
angles  of  attack  at  low  subsonic  speeds. 

Tail  buffeting  received  renewed  interest  during  the  last  decades  with  the  introduction 
of  fighter  aircraft  like  F /  A-18  which  have  high  angle  of  attack  manoeuvring  capability 
typically  in  the  range  of  up  to  35°  but  some  manoeuvres  can  be  performed  at  60°  or 
even  greater.  The  aerodynamics  of  present  day  fighter  aircraft  at  high  incidence  is 
dominated  by  vortex  flows  since  the  wings  are  usually  thin  with  sharp  edges  and  the 
forebody  is  slender.  The  leading  edge  extension  (LEX)  is  used  in  the  design  of  the 
F/A-18  and  other  similar  type  of  aircraft  for  the  purpose  of  generating  additional 
vortex  lift  at  high  angles  of  attack. 

Large  buffet  loads  and  associated  fatigue  damage  have  been  observed  on  the  tails  of 
aircraft  such  as  the  F/A-18  during  certain  flight  conditions,  including  rapid 
maneuvering  at  subsonic  Mach  numbers  and  dynamic  pressures  considerably  less  than 
the  maximum  allowable  values,  see  Scanlon  and  Prey  [2]  and  Streber  and  Rioux  [3]. 
The  cause  of  the  buffet  was  determined  to  be  pressure  fluctuations  within  the  turbulent 
wakes  formed  by  the  breakdown  of  vortices  over  the  upper  surface  of  the  wing  and/or 
leading  edge  extension  (LEX)  surfaces,  see  Lee,  et  al.  [4]  Lee  and  Valerio  [5]  and 
Thompson  [6],  During  high  angle  of  attack  flight  conditions,  the  pressure  fluctuations, 
while  essentially  broad  band  and  chaotic,  revealed  high  levels  of  energy  in  a  narrow 
band  of  frequency.  The  dominant  frequency  corresponding  to  the  spectrum  peak  was 
determined  to  be  very  near  to  the  frequency  of  the  primary  modes  of  tine  tail  structure. 
Consequently  the  structure  was  excited,  producing  larger  than  expected  deflections. 

Attempts  to  alleviate  severe  empennage  buffeting  led  to  a  variety  of  experimental  and 
analytical  investigations  which  attempted  to  clarify  and  quantify  tine  tail  buffet 
phenomenon,  see  Zimmerman  and  Ferman  [7]  and  Zimmerman,  et  al.  [8],  The 
behaviour  of  these  leading  edge  extension  vortices  appeared  to  depend  on  the  angle  of 
attack,  wing  leading  and  trailing  edge  flap  settings,  and  the  presence  of  downstream 
airframe  components  such  as  stabilators  and  vertical  tails  and  to  a  lesser  extent  on  the 
angle  of  sideslip.  Due  to  the  size  of  the  LEXs  and  their  proximity  to  the  forebody,  a 
strong  interaction  can  occur  at  high  angles  of  attack  between  the  LEX  and  forebody 
vortices,  and  small  forebody  contour  distortions  can  promote  large  changes  in  the 
behaviour  of  forebody-LEX  vortex  interactions.  For  example,  small  flow  disturbances 
generated  by  a  nose  boom  can  be  amplified  downstream  by  the  more  powerful  LEX 
vortices  affecting  aircraft  lateral  stability  at  high  angles  of  attack  [9]. 

Separated  flow  from  the  main  wings  has  a  lesser  influence  on  the  vertical  tail,  but  die 
horizontal  stabilators,  being  in  the  wake  of  the  wings,  are  influenced  to  a  much  larger 
extent  than  by  die  forebody  and  LEX  vortices  [4],  The  horizontal  tail  buffet  intensity 
depends  on  the  angle  of  attack  and  stabilator  settings,  and  tiiese  angles  are  related  by 
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Hie  flight  control  schedule.  Detailed  investigation  of  LEX  vortices  on  the  horizontal 
stabilator  has  not  been  carried  out  yet.  However,  based  on  the  F/A-18  IFOSTP  flight 
tests  [3,  10],  excitation  of  the  horizontal  stabilator  was  detected  in  a  range  between  12 
and  28  degrees  of  angle  of  attack. 

hr  order  to  alleviate  the  F/A-18  buffet  problem,  McDonnell  Douglas  Corporation 
focused  their  efforts  on  quantifying  the  buffet  loads  and  response  through  wind  tunnel 
studies  and  flight  tests  [2,  7,  8].  Based  on  some  of  the  studies,  an  interim  solution  for 
the  F/A-18  consisted  of  structural  modifications  to  the  vertical  tail  in  an  attempt  to 
reduce  tire  dynamic  stresses  in  critical  areas.  However,  tire  dynamic  stresses  were  still 
too  severe.  This  led  to  tire  investigation  of  a  LEX  fence  which  is  a  flat  tr  apezoidal  plate 
fixed  perpendicular  to  the  LEX  upper  surface  and  aligned  parallel  to  the  aircraft's 
longitudinal  axis.  The  interaction  between  the  fence  and  the  flow  beneath  the  LEX 
vortex  causes  a  reduction  in  the  dynamic  loads  applied  to  the  tail  and  reduces  the  peak 
accelerations  in  the  first  modes,  but  compromises  the  aircraft's  high  angle  of  attack 
performance  and  maneuverability  [4,  6],  The  development  of  the  LEX  fence  relied 
exclusively  on  hundreds  of  hours  of  costly  wind  tunnel  testing  and  clearly  showed  that 
more  cost  effective  solutions  for  predicting  empennage  buffet  are  required. 

The  importance  of  vertical  tail  buffeting  on  the  fatigue  life  of  high-performance  fighter 
aircraft  was  recognised  by  NASA,  and  tail  buffet  became  one  of  the  research  topics  in 
its  High  Alpha  Technology  Program.  Collective  effort  from  several  US  research  centres 
was  directed  towards  investigation  of  the  behaviour  of  fighter  aircraft  in  the  high  angle 
of  attack  regime  with  the  aim  of  achieving  superior  manoeuvrability  in  future  aircraft 
design.  The  F/A-18  was  chosen  as  a  test  bed  because  it  provides  a  typical  example  of 
the  dominant  vortical  flow-field  (forebody  and  LEX  vortices)  generated  by  future 
advanced  fighter  aircraft  [4].  Tire  program  involved  numerous  flight  tests,  wind  tunnel 
investigations  and  numerical  studies. 

The  Technical  Cooperation  Program  (TTCP),  comprising  the  Defence  Departments  of 
five  nations  (Australia,  Canada,  New  Zealand,  United  Kingdom  and  United  States), 
tasked  one  of  its  technical  panels  (AER  TP5)  to  initiate  a  cooperative  study  on  the 
F /  A-18  at  high  angles  of  attack.  This  study,  that  consisted  of  flight  tests,  wind  tunnel 
and  water  tunnel  tests  of  various  scale  models,  as  well  as  numerical  investigations,  has 
generated  a  large  aerodynamic  database,  that  contains  a  wealth  of  information  on 
aircraft  empennage  buffet  and  buffeting. 

In  1989,  Canada  and  Australia  entered  into  a  collaborative  full-scale  structural  testing 
program,  the  F/A-18  International  Follow  On  Structural  Test  Project  (IFOSTP).  Tire  aim 
of  IFOSTP  was  to  determine  tire  economic  and  safe  life  of  the  F /  A-18  structure  under 
test  loading  which  is  representative  of  Royal  Australian  Air  Force  (RAAF)  and 
Canadian  Forces  (CF)  operation  conditions.  One  of  the  objectives  of  the  IFOSTP  was  to 
develop  methodologies  to  alleviate  vertical  tail  buffet  loads  for  the  F/A-18  life 
extension  program.  A  range  of  wind  and  water  tunnel  experiments  [6,  11,  12]  on 
F/  A-18  sub-scale  models  as  well  as  several  flight  test  programs  [3,  10,  13]  were 


4 


DSTO-RR-0256 


completed  in  support  of  IFOSTP  which  also  provide  valuable  information  for  vertical 
tail  buffeting  studies. 


2.1  Past  Experimental  Studies 

The  appearance  of  empennage  low  speed  buffet  loads  is  known  to  be  caused  by  the 
leading  edge  extension  separation-induced  vortex  breakdown  that  plays  a  major  role  in 
dictating  the  spectral  content  of  the  pressure  fluctuations  on  the  tail  surfaces. 
Characteristics  of  vortex  breakdown  flow  are  complex  and  clearly  nonlinear  in  nature. 
In  the  early  stages  of  an  aircraft  development  process,  the  engineer  often  relies  on 
either  semi-empirical  or  linear  analytic  methods.  These  methods  provide,  at  best,  rough 
estimates  of  frequency  spectrum  distributions  based  on  data  from  similar 
configurations.  Currently  available  prediction  methods  are  inadequate  for  rapid  and 
accurate  prediction  of  specific  time  dependent  buffet  pressure  distributions.  Therefore, 
a  heavy  reliance  on  wind  tunnel  testing  of  small-scale  models  persists. 

Tests  conducted  on  the  F/A-18  have  resulted  in  a  large  database  for  various  model  scales, 
from  small-scale  wind  and  water  tunnels  up  to  full-scale.  Test  data  included  tail  forces 
and  moments,  steady  and  unsteady  pressures,  vortex  breakdown  location,  tip 
accelerations  and  root  strains. 

Flow  visualisation  using  water  tunnels  has  provided  much  sought  after  information  on 
flow  physics  that  is  difficult  to  obtain  from  wind  tunnel  or  flight  tests.  The  results  of 
surface  and  off-body  flow  visualization  information  has  been  found  extremely  useful  in 
the  design  and  implementation  of  control  strategies  to  alleviate  tail  buffet  loads. 


2.1.1  Sub-scale  Model  Testing 

The  largest  amount  of  wind  tunnel  testing  on  the  F/A-18  before  die  TTCP  collaborative 
program  and  IFOSTP  were  established  had  been  carried  out  by  McDonnell  Douglas 
during  the  aircraft  development  stage.  Even  die  development  of  the  LEX  fence 
involved  many  hours  of  wind  tunnel  tests.  Despite  the  substantial  amount  of 
knowledge  on  vertical  tail  buffeting  accumulated  by  McDonnell  Douglas,  most  of  die 
results  are  documented  in  company  reports  and  are  not  available  in  the  open  literature. 
It  is  known  that  the  scale  of  the  McDonnell  Douglas  wind  tunnel  model  is  12%,  and 
rigid  fin  pressures  were  measured  and  scaled  to  full  aircraft  levels  [7,  8],  Some  tests 
were  also  performed  with  a  dexible  tail  tiiat  was  dynamically  scaled  to  match  die 
vibration  frequencies  and  mode  shapes  of  the  full-scale  aircraft.  Strain  gauges  were 
installed  on  the  wind  tunnel  model  for  bending  and  torsion  moment  measurements, 
and  accelerometers  on  the  fin  tip  to  detect  the  response  levels.  The  scaling  was  carried 
out  such  that  the  dynamic  pressures  for  the  wind  tunnel  tests  simulate  those  at  Eight. 

Some  of  the  first  experimental  investigations  of  die  vertical  tail  buffet  of  the  F/A-18 
models  have  been  conducted  by  Sellers,  at  al.  [14],  Erickson,  at  al.  [15],  Wentz  [16]  and 
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Lee  and  Brown  [17].  These  experiments  showed  that  die  vortex  produced  by  the  wing 
LEX  breaks  down  ahead  of  the  vertical  tails  at  angles  of  attack  of  25°  and  higher  and 
the  breakdown  flow  produced  unsteady  loads  on  the  vertical  tails. 

Erickson  [15]  documented  results  of  a  wind  tunnel  test  of  a  6%-scale  model  of  the 
F/A-18.  He  concluded  that  the  high  Reynolds  number  behaviour  of  the  LEX  vortex 
flow  can  be  simulated  at  lower  Reynolds  numbers  in  the  wind  tunnel  and  diat  the 
vortex  core  breakdown  is  relatively  unaffected  by  Mach  number  until  shock  waves 
appear  over  the  wings  that  interact  with  the  vortices.  He  also  found  that  the  F/  A-18 
forebody  vortices  are  comparatively  weak  and  are  dominated  by  the  LEX  vortices  at 
all  Mach  numbers. 

Wentz  [16]  performed  water  tunnel  tests  to  study  the  flow  associated  with  the  tail  buffet  of  the 
F/A-18  aircraft  using  a  l/48dl  scale  model.  Coloured  dyes  were  used  for  flow  visualization  of 
LEX  vortex  breakdown  to  determine  vortex  patterns.  He  also  used  surface  hot  film  gauge 
placed  on  the  outboard  surface  of  the  starboard  tail  to  measure  the  heat  transfer  at  the 
surface  and  to  study  the  turbulent  energy  and  the  frequencies  present  in  the  unsteady 
wake  flow.  Hot  film  anemometer  signals  showed  that  the  tail  surface  turbulence  increased 
with  angle  of  attack,  and  that  dominant  frequencies  appear  in  the  flow  when  breakdown 
occurred.  Comparisons  with  available  wind  tunnel  measurements  showed 
approximately  similar  values  of  the  vortex  trajectories,  breakdown  positions  and 
Strouhal  numbers  for  unsteady  fin  pressures.  These  confirm  the  assumption  that  the 
fundamental  flow  patterns  are  independent  of  Reynolds  number  and  the  water  tunnel 
results  are  reliable  when  used  to  supplement  wind  tunnel  and  flight  test  data. 

Investigations  of  vertical  tail  buffeting  and  the  effect  of  the  LEX  fence  on  buffet  loads 
were  performed  by  Lee  and  Brown  [17]  on  a  rigid  6%  scale  model  of  the  F/A-18.  The 
starboard  vertical  tail  was  instrumented  with  24  fast  response  absolute  pressure 
transducers  on  each  of  the  inboard  and  outboard  surfaces.  Several  strain  gauges  and 
accelerometer  were  installed  at  the  starboard  tail  and  tire  port  tail  was  also 
instrumented  with  strain  gauges  to  provide  accurate  measurements  of  the  normal 
force,  root  bending  and  torsion  moments.  The  studies  also  included  the  forebody 
pressure  distributions,  steady  and  unsteady  pressure  measurements  on  the  LEX  and 
wing,  as  well  as  the  flowfield  behind  the  vertical  tails.  Spectral  analyses  of  the  pressure 
fluctuations  on  die  tail  mid  in  the  vortical  flow  region  behind  die  tail  revealed  the 
presence  of  dominant  oscillations  during  buffet.  Large  reductions  in  steady  mid 
unsteady  pressures  were  detected  on  die  vertical  tail  inboard  surfaces  at  30  degrees 
angle  of  attack  witii  the  LEX  fence  installed. 

Hebbar,  Platzer  mid  Cavazos  [18]  studied  the  effect  of  pitch  rate  on  the  development 
of  vortex  breakdown  using  flow  visualization  results  of  a  2%  F/A-18  model  tested  in 
a  water  tunnel.  Results  indicated  the  presence  of  hysteresis  in  the  instantaneous 
position  of  the  vortex  breakdown  point  relative  to  the  static  case.  Thus,  the  vortex 
breakdown  point  moves  rearward  with  increased  pitch-up  motion  and  forward  with 
increased  pitch-down  motion.  Vortex  core  breakdown  position  was  also  determined 
to  be  a  function  of  yaw  angle.  The  windward  vortex  breakdown  location  was  seen  to 
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move  forward  and  inboard,  while  the  leeward  vortex  breakdown  point  travelled  aft 
and  outboard.  The  results  of  the  study  confirmed  that  aircraft  flight  dynamics  and 
vortex  breakdown  dynamics  must  be  considered  as  coupled  events. 

Lee  and  Tang  [19]  measured  the  tail  buffet  on  a  rigid  6%-scale  wind  tunnel  model  of 
the  F/  A-18.  Strain  gauges  were  used  to  measure  force  and  moments  on  the  vertical 
tail  surface.  Fast  response  pressure  transducers  were  also  installed  in  the  tail  surface. 
Comparison  of  steady  forces  and  moments  on  the  vertical  tail  determined  from 
pressure  summation  and  strain  gauge  measurements  showed  good  agreement.  They 
concluded,  that  probability  density  functions  of  the  vertical  tail  loads  as  well  as 
pressures  on  the  tail  surface  at  low-pitch  frequencies  could  be  represented  by 
Gaussian  distributions. 

An  experimental  investigation  of  vortex-tail  interaction  on  a  76°  sharp-edged  delta 
wing,  twin  vertical- tail  configuration  has  been  conducted  by  Washburn,  et  nl.  [20].  The 
vertical  tails  were  placed  at  nine  locations  behind  the  delta  wing.  The  experimental 
data  showed  that  the  aerodynamic  loads  are  more  sensitive  to  the  chordwise  tail 
location  than  its  spanwise  location.  Although  the  tail  location  did  not  affect  the  vortex 
core  trajectories,  it  affected  the  location  of  vortex-core  breakdown,  the  aerodynamic 
loads,  and  the  tail  buffeting  levels.  As  the  tails  were  moved  laterally  toward  the  vortex 
core  trajectory,  the  buffeting  response  and  excitation  were  reduced.  Moreover,  the 
investigation  showed  that  the  presence  of  a  flexible  tail  could  affect  the  unsteady 
pressures  on  the  rigid  tail  on  the  opposite  side  of  the  model.  Although  the 
configuration  was  less  complex  than  the  F/  A-18,  results  showed  similarities  of  buffet 
flow  characteristics.  As  such,  this  case  presents  an  opportunity  to  isolate  the  primary 
flow  features  of  concern  to  study  the  phenomenon  and  provides  data  for  validation 
of  analytical  methods  using  a  simplified  model. 

A  buffet  research  program  was  initiated  at  NASA  Langley  Research  Center  (LaRC)  to 
study  vortex  flowfield  mechanisms,  with  a  focus  on  interaction  of  the  vortical 
flowfields  on  tail  buffet  stability  and  control,  and  develop  a  database  relating 
configuration  features  to  these  characteristics.  Topics  that  were  directly  related  to  tail 
buffet  modelling  and  prediction  included  pressure  and  response  measurements,  and 
the  effects  of  variation  in  geometry  on  vortex  flow  and  subsequent  impact  on  tail  buffet 
loads.  Shah,  et  al.  [21]  reported  on  to  effects  of  leading  edge  extension  modifications  on 
the  F/A-18  vertical  tail  buffet  loads.  He  used  a  16%  scale  model  of  the  production 
F /  A-18  designed  originally  as  a  free  flight  drop  model.  The  port  vertical  tail  was  rigid 
while  the  starboard  tail  was  flexible  and  dynamically  scaled  to  match  the  stiffness, 
inertia  and  natural  structural  frequencies  of  the  F/A-18  full-scale  aircraft  vertical  fins. 
The  buffeting  characteristics  on  the  flexible  tail  were  acquired  in  terms  of  the  tail  root 
bending,  mid-span  torsion,  and  leading  edge  acceleration.  Some  flow  visualisation  was 
conducted  using  a  laser  light  sheet  technique. 

Moss,  Cole  and  Doggett  [22]  also  tested  a  rigid  16%,  full-span  model  of  the  F/A-18 
aircraft  that  was  fitted  with  flexible  vertical  tails  of  different  stiffness.  Two  flexible  tails 
were  constructed  which  had  the  same  planform  geometry  as  the  full-scale  F/A-18,  but 
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were  not  dynamically  scaled.  Vertical-tail  buffet  response  results  were  obtained  over 
the  range  of  angle  of  attack  from  -10°  to  +40°,  and  over  the  range  of  Mach  numbers 
from  0.3  to  0.95.  The  results  were  presented  in  the  form  of  the  buffet  excitation 
parameter  and  indicated  that  the  buffet  response  occurs  in  the  first  bending  mode, 
increases  with  increasing  dynamic  pressure  and  is  larger  at  M=0.3  than  that  at  a  higher 
Mach  number. 

Moses,  et  al.  [23,  24]  reported  results  for  a  16%  scale  wind  tunnel  model  of  the  F/A-18 
employed  to  study  spatial  correlation  of  the  buffeting  flow  at  high  angles  of  attack. 
Unsteady  pressures  were  measured  on  rigid  and  flexible  tails.  Cross-correlation  and 
cross-spectral  density  functions  were  presented  which  illustrate  the  time  lags  and  phase 
lags  associated  with  the  unsteady  differential  pressures  at  stations  on  die  tails.  The  time 
and  phase  lags  are  characteristic  of  a  wave  and  were  shown  to  be  functions  of  die 
distance  between  stations  and  the  transport  velocity.  At  a  given  angle  of  attack,  the 
correlation  reportedly  scaled  widi  speed,  as  determined  through  comparisons  with  data 
from  odier  documented  tests.  Tail  dexibility  did  not  appear  to  significantiy  affect  the  time 
or  phase  delays.  The  results  presented  by  Moses  provide  a  better  understanding  of  the 
complex  phenomenon  of  buffeting,  and  the  statistical  analysis  showing  the  existence  of 
coherent  structures  in  the  flow  over  the  tail's  surface  is  useful  in  devising  active  buffet 
load  alleviation  technologies. 

The  initial  dow  visualisation  tests  performed  at  the  Platform  Sciences  Laboratory  of  the 
Defence  Science  and  Technology  Organisation  (DSTO)  were  devoted  mainly  to 
measurements  of  the  vortex  burst  positions  and  comparisons  with  water  tunnel  and 
results  from  other  wind  tunnel  studies,  see  [11].  It  was  shown  by  Thompson  [11]  that 
die  vortex  burst  position  was  not  affected  by  the  dap  settings,  horizontal  stabilator 
settings,  model  size  and  LEX  fence  installation.  However,  he  found  that  the  presence  of 
engine  inlet  dow  affects  the  behaviour  of  die  LEX  vortex  and  tends  to  move  the  burst 
point  downstream  for  a  fixed  angle  of  attack.  Refinement  of  dow  visualisation 
techniques  reported  by  Thompson  [6]  allowed  studying  the  topology  of  die  streamlines 
in  die  vicinity  of  die  LEX  fence.  This  study  shed  some  light  on  the  effect  of  the  LEX 
fence  on  die  vortical  dow  field. 

Martin,  et  al.  [12]  presented  results  of  buffet  measurements  carried  out  on  a  1/9*  scale 
model  of  die  F/A-18  in  9ft  x  7ft  low-speed  wind  tunnel  at  die  DSTO.  The  model  was 
mounted  on  a  pitch/ roll  rig  via  a  six-component  strain  gage  balance  and  has 
dow-through  engine  inlet  simulation.  A  fast  response  pressure  transducer  and  an 
accelerometer  were  located  on  die  outboard  side  of  the  port  fin.  To  measure  differential 
steady  pressures,  there  were  tiiree  pairs  of  pressure  taps,  located  on  opposite  sides  of 
die  fin.  A  ground  vibration  test  was  performed  on  die  model  fin  to  determine  die  fin 
primary  bending  mode  shapes  and  the  natural  frequencies.  The  first  natural  bending 
was  found  to  be  68  Hz,  which  does  not  correspond  to  the  aircraft  frequency  of  15  Hz 
since  die  model  was  not  aeroelastically  scaled.  However,  excitation  of  die  fin  would  be 
expected  since  diis  value  fell  widiin  the  range  of  frequencies  that  occur  in  the  vortex 
burst  of  die  1/9*  scale  model  at  the  tunnel  speeds  where  the  tests  were  conducted. 
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Results  of  Martin,  et  al.  provided  valuable  insight  into  the  character  of  the  vortex 
breakdown  and  subsequent  empennage  vibration  on  the  F/A-18.  Flow  visualization 
measurements  with  the  vertical  tail  on  and  off  showed  that  the  presence  of  the 
vertical  tail  moves  the  point  of  breakdown  upstream.  Pressure  measurements  on  the 
wing  surface  beneath  the  vortex  breakdown  and  on  the  vertical  tail  showed  that 
vortex-breakdown  pressure  field  contained  energy  over  a  moderately  narrow 
frequency  band.  The  characteristic  frequency  was  seen  to  be  a  linear  function  of  free 
stream  velocity  and  to  vary  with  angle  of  attack. 


2.1.2  Full  Scale  Tail  Buffet  Tests 

Full-scale  tail  buffet  tests  at  the  NASA  Ames  Research  Center  were  performed  in  order 
to  study  die  flowfield  characteristics  of  tail  buffet  over  a  wide  range  of  angles  of  attack 
and  sideslip,  to  quantify  the  effects  of  LEX  fence  in  reducing  the  tail  buffet  loads,  and 
to  provide  full-scale  data  for  comparison  with  sub-scale  wind  tunnel  model  results. 
The  test  article  was  a  full-scale  F/A-18  supplied  by  the  US  Navy  with  the  aircraft 
engines  and  avionics  removed.  The  tests  were  conducted  in  the  80  x  120  ft  wind  tunnel 
of  the  National  Full-Scale  Aerodynamic  Complex.  The  maximum  dynamic  pressure 
attainable  was  33psf  at  a  maximum  velocity  of  100  knots. 

In  tire  first  entry,  only  32  transducers  were  installed  on  the  port  vertical  tail  in  a  4  x  4 
grid  on  each  side  of  the  tail.  Four  accelerometers  were  installed  near  the  tail  tips  to 
record  response  in  botii  bending  and  torsion  modes.  The  angle  of  attack  ranged  from 
18  to  50  degrees  and  the  sideslip  varied  from  -15  to  15  degrees.  Results  on  bending 
moment  variations  with  angle  of  attack,  individual  and  differential  pressure  power 
spectral  densities,  RMS  pressure  coefficients  and  peak  power  frequencies  for  the  first 
test  were  reported  by  Meyn,  et  nl.  [25],  Meyn  and  James  [26, 27]. 

In  the  second  wind  tunnel  entry,  pressures  were  measured  using  a  finer  6x8  grid  on 
each  side  of  the  tail  and  James  and  Meyn  [28]  showed  the  effect  of  spatial  resolution  of 
pressure  measurements  on  integrated  buffet  loads.  Their  results  indicate  that  the 
substantial  difference  in  the  results  between  the  two  wind  tunnel  entries  was  due  to 
insufficient  number  of  transducers  used  in  the  first  test. 

Another  entry  in  tire  80  x  120  ft  wind  tunnel  was  conducted  by  the  USAF  [29],  For  tire 
majority  of  the  test  runs,  the  conditions  were  similar  to  those  used  in  the  Meyn  and 
James  [28]  studies.  However,  only  72  transducers  were  installed  on  the  starboard 
vertical  tail.  A  limited  number  of  runs  were  carried  out  at  a  dynamic  pressure  of  20  psf 
with  and  without  the  LEX  fence.  The  angle  of  attack  was  varied  from  20  to  40  degrees 
and  the  sideslip  from  -16  to  16  degrees.  Comparison  of  the  RMS  differential  pressure 
coefficient  maps  showed  that  the  LEX  fence  not  only  reduces  the  overall  magnitude  of 
the  buffet  fluctuations,  but  also  alters  the  shape  of  tine  pressure  contours  along  tine 
leading  edge  and  in  the  tip  region. 
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Meyn,  et  al.  [30]  presented  a  correlation  of  power  spectral  characteristics  of 
differential  pressures  measured  in  flight  and  on  a  full-scale,  production  F/A-18 
aircraft  in  the  wind  tunnel.  For  the  full-scale  F/A-18  wind  tunnel  test,  the  maximum 
wind  speed  obtained  created  a  spectrum  of  differential  pressures  well  below  the  first 
bending  mode  of  the  vertical  tail.  Therefore,  the  tails  were  considered  rigid. 
Comparisons  of  power  spectra  of  the  rigid  model  pressure  coefficients  with  power 
spectra  of  flight  data  indicated  that  tail  flexibility  resulted  in  lower  pressure 
coefficient  values.  Further  comparisons  were  made  with  results  from  rigid  tails  on  a 
12%-scale  wind  tunnel  model.  For  angles  of  attack  between  28  and  36  degrees,  the 
peak  power  pressure  coefficient  was  highest  on  the  12%  model  and  lowest  on  the 
aircraft.  This  trend  seems  to  indicate  that  the  highest  pressures  are  obtained  on  tails 
that  are  most  rigid. 

Results  of  full-scale  tail  buffet  investigations  at  NASA  Ames  were  also  compared  with 
those  of  wind  tunnel  tests  on  a  16%-scale  model  and  Moses  and  Pendleton  [31]  found 
drat  the  results  from  the  two  models  revealed  similarities  in  the  trends  of  the  spectral 
contents  as  a  function  of  angle  of  attack.  They  concluded  that  characteristics  of  the 
differential  pressures  on  tire  fin  surface  on  smaller  wind  tunnel  models  could  be  scaled 
up  to  full  scale  with  good  accuracy. 


2.2  Numerical  Investigation  of  Tail  Buffeting 

As  future  military  aircraft  are  also  expected  to  operate  at  high  angles  of  attack  in 
maneuvering  flight,  practical  methods  are  needed  to  accurately  account  for  the  buffet 
loads  early  in  air  aircraft  development  process.  Present  day  methods  rely  mainly  on 
semi-empirical  analysis  and  use  of  linear  aerodynamics  to  establish  the  necessary 
buffet  loads  and  tail  dynamic  response  for  high  angle  of  attack  conditions.  A  popular 
method  developed  by  Zimmerman  and  Ferman  [7,  8]  is  based  on  tire  premise  that  some 
form  of  buffet  data  is  available  through  sub-scale  wind  tunnel  tests. 

However,  tire  use  of  linear  aerodynamics  is  not  quite  correct  since  the  results  from 
flight  tests  on  tire  F/A-18  vertical  tails  show  large  regions  of  separated  flow.  Tire  idea 
of  using  linear  aerodynamics  is  based  on  the  assumption  that  even  though  the  flow  on 
the  rigid  model  is  highly  nonlinear  with  regions  of  flow  separation,  the  perturbations 
on  the  flow  due  to  the  tail  vibrations  can  still  be  approximately  determined  from  linear 
theory.  This  may  be  valid  for  attached  flows  with  small  amplitudes  of  tail  motion,  but 
when  flow  separation  is  present,  the  tail  motion  will  certainly  change  tire  separation 
and  attachment  lines  on  the  tail  surfaces.  The  study  of  aeroelasticity  in  separated  flows 
is  a  challenging  problem  and  Forsching  [32]  discussed  tire  use  of  semi-empirical 
aerodynamics  for  separated  flows  in  estimating  the  structural  response  during 
buffeting. 

Much  work  is  required  before  a  better  aerodynamic  model  can  be  substituted  for  tire 
widely  used  linear  theory.  There  is  presently  no  analytical  method  for  the  prediction  of 
unsteady  separated  flows  during  buffeting.  Tire  use  of  computational  fluid  dynamics 
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(CFD)  is  promising  but  the  computing  cost  is  still  prohibitively  expensive.  The  only 
alternative  is  to  test  dynamically  scaled  flexible  models  in  the  wind  tunnel,  and  that  is 
also  an  extremely  expensive  and  time  consuming  process.  Ideally,  this  requires 
matching  the  Reynolds  number  and  Mach  number  exactly  to  satisfy  scaling  laws  [7,  8] 
in  order  that  measurands  from  a  wind  tunnel  environment  can  be  converted  to  actual 
flight  conditions. 


2.2.1  Linear  Aeroelastic  Models 

As  the  buffet  problem  is  inherently  aeroelastic  in  nature,  its  solution  depends  on  the 
advances  in  theoretical  research  in  the  areas  of  structural  dynamics  and  unsteady 
aerodynamics  as  well  as  numerical  solution  schemes  associated  with  computational 
aeroelasticity.  General  aeroelastic  phenomena  were  described  by  Theodorsen  [33], 
Bisplinghoff,  Adhley  and  Halfman  [34],  and  Dowell  [35],  A  more  complete  formulation 
of  the  general  aeroelastic  equations  was  described  by  Bisplinghoff  [36]  and  Fung  [3 7], 
Pines  [38]  presented  physical  models  of  flutter,  examined  the  governing  equations  in  a 
general  sense  and  established  relationships  for  which  flutter  can  exist. 

In  the  1960s,  only  simple  computational  aeroelastic  models,  such  as  a  two-dimensional 
wing  mounted  on  an  elastic  support  or  lumped  mass  structural  models  with 
theoretical,  or  semi-empirical  aerodynamics,  were  applied  in  aeroelastic  design.  The 
first  aeroelastic  analyses  considered  incompressible  flow  and  linear  springs.  Although 
convenient  for  their  simplicity  and  minimal  computational  requirements,  they  lack  the 
reality  of  three-dimensional  nonlinear  flows.  Beginning  in  the  1970s,  structural  analysis 
codes  based  on  the  Finite  Element  Method  (FEM)  initiated  the  use  of  aero-structural 
models  by  implementing  flat  plate  vortex  lattice  methods.  While  this  provided  for 
three-dimensional  effects,  the  analysis  was  restricted  to  linear  aerodynamics,  often 
neglecting  thickness  effects. 

In  the  early  1980s,  a  significant  advance  was  made  when  the  analysis  was  extended  to 
transonic  flows.  Because  transonic  flows  are  described  by  nonlinear  equations,  there 
were  methods  developed  for  simultaneous  numerical  integration  of  the  governing 
equations  for  the  flow  and  the  dynamic  equations  of  motion  of  tire  structure.  The  first 
efforts  in  this  pioneering  work  were  by  Rizetta  [39,  40],  Ballhaus  and  Goorjian  [41], 
Yang,  Guruswamy,  and  Striz  [42],  Yang,  et  al.  [43],  and  Guruswamy  and  Yang  [44], 
Their  work  was  the  first  to  truly  capture  dynamic/ aerodynamic  interaction  by 
considering  the  flowing  fluid  and  airfoil  together  to  be  a  single  dynamic  system. 

Desmarais  and  Bennett  [45]  described  an  approach  for  flutter  analysis  where  tire 
generalised  mass  and  mode  shapes  are  determined  through  the  use  of  an  external 
structural  dynamics  analysis.  These  modes  describe  the  shape  of  the  wing  and  are  used 
to  compute  an  aerodynamic  influence  matrix.  The  aerodynamics  are  modelled  using  a 
subsonic  kernel  function  approach.  The  solution  is  performed  in  the  frequency  domain, 
and  the  roots  of  the  eigensolution,  which  are  dependant  upon  tine  velocity,  density,  and 
Mach  number,  are  used  to  determine  the  critical  conditions  necessary  for  flutter.  The 
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governing  equations  are  cast  in  a  linear  form,  and  a  simple  harmonic  form  of  the 
solution  is  assumed  and  application  of  the  method  is  limited  by  these  assumptions. 

Frequency  domain  solutions  are  the  most  widely  used  since  they  require  less  resources 
and  are  a  more  convenient  solution  scheme.  However,  a  frequency  domain  solution 
does  not  provide  a  physical  description  of  the  motion  as  do  time  domain  solutions 
which  provide  a  more  general  approach,  and  do  not  assume  a  form  of  the  solution,  i.e. 
harmonic  motion,  as  do  the  frequency  domain  schemes. 

One  of  the  first  time  domain  solutions  of  the  governing  aeroelastic  equations  was 
developed  by  Devers  [46].  He  used  a  vortex-lattice  method  as  an  aerodynamic  model. 
However,  simplifications  in  his  model  prevented  the  unsteady  nature  of  trailing  edge 
and  wing-tip  vortex  effects  from  being  considered.  In  addition,  his  integration  scheme 
assumed  a  form  of  the  solution  and,  consequently,  the  simulation  was  limited  by  this 
assumption. 

In  the  early  1980s,  Guruswamy  and  Goorjian  [47]  successfully  demonstrated  an 
unsteady  transonic  small  disturbance  CFD  code  that  was  combined  with  a  modal 
analysis  computational  structural  dynamics  code  to  calculate  the  transonic  flutter  of  an 
unswept  rectangular  wing  which  was  later  extended  to  include  the  effects  of  moving 
control  surfaces. 

Batina  [48,  49]  developed  a  method  that  solved  the  inviscid  transonic  small  disturbance 
(TSD)  equations  for  use  in  aeroelastic  analysis.  This  work  was  undertaken  in  order  to 
provide  an  alternative  to  the  widely  used  linear  methods,  such  as  doublet  lattice,  and 
to  improve  upon  previous  TSD  codes.  The  code  developed  to  solve  the  TSD  equation 
allowed  transonic  effects  to  be  determined  in  flutter  calculations.  The  method  was 
limited  to  inviscid,  irrotational,  attached  flows  and  is  invalid  for  the  prediction  of  high 
angle  of  attack  vortex-induced  tail  buffet. 

Some  attention  has  also  been  given  to  the  problem  of  modelling  the  aeroelastic 
response  of  flexible  wings.  Among  the  earliest  investigators  to  capture 
dynamic /aerodynamic  interaction  by  considering  the  wing  and  the  flowing  fluid  to  be 
a  single  dynamic  system  were  Eastep  and  Olsen  [50]  and  Borland  and  Rizzetta  [51]. 
These  efforts  were  directed  toward  transonic  flows  and  used  small  disturbance  models 
of  the  flowfield. 


2.2.2  CFD/ CSD  Aeroelastic  Solvers 

By  the  1990s  advances  in  aeroelastic  CFD  and  computer  speeds  led  to  a  more 
sophisticated  analysis  tool  through  the  introduction  of  three-dimensional  unsteady 
CFD  solvers  combined  with  a  computational  structural  dynamics  (CSD)  modal  model 
analysis. 
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Strganac  [52]  and  Strganac  and  Mook  [53]  also  analysed  the  unsteady  aeroelastic 
behaviour  of  flexible  wings.  They  developed  a  method  for  analysing  the  aeroelastic 
response  of  rigid  and  flexible  finite-aspect-ratio  wings  mounted  on  elastic  supports. 
The  method  is  based  on  the  general  unsteady  vortex-lattice  method  and  can 
accommodate  non-linear  structural  models  and  small,  but  finite,  amplitudes  of  motion 
that  make  it  suitable  for  limit  cycle,  i.e.  post  flutter,  analyses. 

In  their  approach,  Strganac  and  Mook  followed  a  procedure  that  was  first  proposed  by 
Loring  [54],  who  introduced  generalised  coordinates  through  the  use  of  the  natural 
modes  to  solve  the  flutter  problem.  They  expressed  the  deflections  as  expansions  in 
terms  of  the  linear  free-vibration  modes  and  considered  the  time-dependant 
coefficients  in  the  expansions  to  be  the  generalised  coordinates  of  the  motion.  The 
aerodynamic  loads  were  also  expressed  as  expansions  in  the  modes.  Following 
Galerkin's  procedure,  they  reduced  the  original  set  of  governing  partial  differential 
equations  to  a  finite  set  of  ordinary  differential  equations.  The  governing  equations  are 
developed  for  motion  about  arbitrary  static  angles  of  attack,  providing  a  capability  to 
study  the  associated  nonlinear  effects. 

The  three-dimensional  unsteady  vortex-lattice  method,  as  used  by  Strganac  and  Mook, 
is  able  to  treat  arbitrary  manoeuvres  of  wings  of  arbitrary  planforms,  including  highly 
swept  delta  wings  which  exhibit  leading-edge  separation,  configurations  of  multiple 
closely  coupled  lifting  surfaces,  and  low  aspect  ratio  planforms.  It  can  also  treat 
arbitrary  angles  of  attack  and  camber,  as  long  as  stall  or  vortex  bursting  in  the  near 
wake  does  not  occur.  The  method  accurately  predicts  spanwise  distributions  of  static 
loads  and  accounts  for  the  nonlinear  effects  of  the  wakes  adjoining  the  tips  and  trailing 
edge.  Because  the  governing  equations  are  solved  in  the  time  domain,  the  method 
provides  an  opportunity  to  analyse  the  results  using  visualisation  of  the  flowfield  and 
wing  motion. 

Schuster,  et  al.  [55]  presented  an  aeroelastic  approach  based  on  a  blocked  structured 
Euler/Navier-Stokes  CFD  method.  Since  that  time,  a  number  of  static  and  dynamic 
rigid  and  aeroelastic  test  cases  have  been  analyzed.  The  Navier-Stokes  solver  was 
coupled  with  a  linear  structural  modal  model  to  provide  static  aeroelastic  analysis  of 
wings  at  extreme  flight  conditions.  The  method  compared  well  with  experimental  data 
and  was  later  extended  to  enable  analysis  of  more  general  configurations  including 
vertical  tail  surfaces. 

Guruswamy  [56]  reported  on  an  Euler  CFD  based  solver  coupled  with  a  modal 
structural  model  for  computing  aeroelastic  responses  of  wings.  Later  he  extended  die 
method  to  include  a  Navier-Stokes  CFD  based  dynamic  aeroelastic  analysis  capability. 
The  computed  results  compared  well  with  the  available  experimental  data  for  steady 
and  unsteady  flows.  Effects  of  flexibility,  wing  sweep  angle  and  pitch  rate  were 
included  in  the  prediction  of  aeroelastic  response  of  a  wing-alone  configuration. 

Kandil,  Kandil  and  Massey  [57]  presented  the  first  computational  simulation  of  the 
vertical  tail  buffet  using  a  delta  wing-vertical  tail  configuration.  A  76°  swept 
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sharp-edged  delta  wing  has  been  used  along  with  a  single  rectangular  vertical  tail 
which  was  placed  aft  of  the  wing  along  tire  plane  of  geometric  symmetry.  The  tail  was 
allowed  to  oscillate  in  bending  modes.  The  flow  conditions  and  wing  angle  of  attack 
have  been  selected  to  produce  an  unsteady  vortex-breakdown  flow.  Unsteady  vortex 
breakdown  of  leading-edge  vortex  cores  was  captured,  and  unsteady  pressure  forces 
were  obtained  on  the  tail.  Their  computational  results  are  in  good  qualitative 
agreement  with  the  experimental  data  of  Washburn,  Jenkins  and  Ferman  [20], 

Kandil,  Massey  and  Sheta  [58]  studied  the  effects  of  coupling  and  uncoupling  the 
bending  and  torsional  modes  and  tire  flow  Reynolds  number  on  the  buffet  response  of 
a  single  rectangular  tail.  It  has  been  shown  that  tire  coupled  response  produces  higher 
deflection  than  that  of  tire  uncoupled  response.  Moreover,  the  response  of  the  coupled 
case  reaches  periodicity  faster  than  drat  of  the  uncoupled  case.  It  has  also  been  shown 
drat  dre  deflections  of  the  low-Reynolds  number  case  are  substantially  lower  than  those 
of  the  high  Reynolds  number  case. 

Kandil,  Sheta  and  Massey  [59]  presented  a  finite  volume  CFD-based  study  including 
buffet  response  of  tire  F/A-18  twin  tails  in  a  generic  delta  wing-fin  configuration.  The 
leading-edge  vortex  breakdown  How  has  been  generated  using  a  76°-swept  back 
sharp-edged  delta  wing  which  was  pitched  at  32°  angle  of  attack.  Each  of  tire  twin  tails 
was  cantilevered  at  a  trailing  edge  extension  of  the  delta  wing.  A  multi-block  grid  was 
used  to  solve  tire  problem  for  two  lateral  locations  of  the  twin  tails.  Only  uncoupled 
beirding-torsion  response  cases  were  considered  in  this  study.  While  validation  against 
experimental  data  was  not  reported,  the  results  represent  expected  trends  showing 
promise  for  tire  metirod. 

Computational  simulations  on  tire  complete  F/A-18  geometry  for  long  time  histories 
required  for  statistical  analysis  are  prohibitively  expensive  and  have  been  carried  out 
only  in  large  research  organisations  such  as  NASA  Ames  and  NASA  LaRC  where 
sufficient  computing  resources  are  available.  One  of  the  first  attempts  at  numerical 
simulation  of  die  F/A-18  tail  buffeting  found  in  die  open  literature  was  undertaken  by 
Rizk,  et  nl.  [60]  and  Gee,  et  nl.  [61].  They  carried  out  unsteady  tiow  calculations  around 
die  F/A-18  aircraft  and  also  investigated  buffeting  of  the  vertical  tails  at  M=0.243  and 
30.3  degrees  angle  of  attack  and  compared  predictions  with  wind  tunnel  tests  to  show 
die  accuracy  of  tiieir  numerical  schemes. 

In  Rizk,  et  al.  [60],  die  computations  were  carried  out  for  a  rigid  vertical  tail,  so  no 
aerodynamic  coupling  was  taken  into  account.  Also,  only  a  symmetrical  configuration 
was  considered.  Comparisons  with  full-scale  F/  A-18  wind  tunnel  tests  showed  rather 
poor  correlation  and  Gee,  et  nl.  [61]  improved  the  CFD  results  using  an  enhanced  grid 
system.  The  new  grid  system  contained  a  finer  grid  resolution  in  die  fuselage  forebody 
and  LEX  regions,  and  included  die  details  of  the  engine  inlet  geometry.  This  led  to 
more  accurate  prediction  of  die  LEX  vortex  burst  position  and  improved  prediction  of 
die  vertical  fin  unsteady  loads. 
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When  compared  with  experiments,  the  predictions  of  Gee,  et  al.  [61]  were  found  to  be 
quite  good  in  some  of  the  pressure  results,  but  were  poor  in  the  fin  root  bending 
moment  and  peak  power  coefficients.  Inclusion  of  the  aeroelastic  deformation  in  the 
numerical  simulations  is  necessary  to  improve  the  predictions,  as  the 
structural-induced  aerodynamic  loads  can  be  significant  at  higher  dynamic  pressures. 
This  requires  simultaneous  solutions  of  the  aerodynamic  loads  and  structural  response 
of  the  vertical  tail  so  that  the  interaction  at  each  time  step  can  be  accounted  for. 

The  most  rigorous  attempt  to  computationally  simulate  the  F/A-18  vertical  fin 
buffeting  was  undertaken  by  Findlay  [62],  He  used  a  tightly  coupled  dynamic 
aeroelastic  analysis  to  predict  empennage  buffet  at  30  degrees  angle  of  attack.  The 
method  employs  the  solution  of  the  time-dependent  Navier-Stokes  equations  to 
describe  the  unsteady  flowfield  and  a  structural  dynamics  modal  analysis  to  compute 
response  of  the  structure.  The  model  was  employed  to  perform  a  computational  study 
of  steady  and  unsteady  calculations  of  the  F/A-18  aircraft  with  rigid  and  flexible  tails. 

The  results  of  Findlay  showed  reasonable  agreement  with  in-flight,  and  full-scale  wind 
tunnel  test  measured  trends.  However,  the  absolute  magnitudes  of  buffet  loads  and 
dominant  frequencies  of  the  buffet  pressures  were  under-predicted.  This  is  primarily 
attributed  to  insufficient  grid  resolution,  as  increases  in  flowfield  field  fidelity  showed 
improvements  in  surface  pressures  for  steady-state  calculations.  However,  as  reported 
in  [62],  computational  requirements,  expected  to  be  adequate  for  unsteady  calculations 
with  a  grid  of  sufficient  fidelity  were  beyond  the  resources  available. 

As  reported  by  Findlay,  the  effect  of  flexibility  of  the  tail  showed  a  decrease  in  buffet 
pressure  magnitude  from  the  values  for  the  rigid  calculations  which  is  consistent  with 
experimental  results.  The  difference  in  the  two  values  compared  well  with  measured 
differences  between  rigid  and  flexible  model  tests.  This  is  an  important  finding  as 
buffet  is  often  assumed  to  be  independent  of  tail  response  when  computing  unsteady 
aerodynamic  loads.  Thus,  rigid  computations  usually  over-predict  tire  flight  data  and 
comparison  between  rigid  computations  and  flexible  flight  test  data  could  be 
misleading. 


3.  Structural  Dynamic  Model 

A  general  approach  to  the  computation  of  the  response  of  a  structure  to  dynamic  loads 
is  provided  by  energy  methods,  see  Meirovitch  [63].  In  such  methods,  the  solution  is 
derived  by  describing  the  dynamic  response  in  terms  of  small  displacements  about 
equilibrium  points. 
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3.1  Governing  Equations  of  Motion 

The  equation  governing  the  motion  of  the  dynamical  system  representing  the  vertical 
tail  can  be  obtained  by  Hamilton's  principle  [63],  which  can  be  stated  in  the  form 

I  L 

J \{dEK  -  BE,,  +  dWA  )dldt  =  0  (1) 

0  0 


where  EK  and  EP  are  kinetic  and  potential  energy  of  the  system,  W A  is  the  work 
performed  by  the  aerodynamic  loading  and  L  is  the  vertical  tail  length. 

To  facilitate  the  development  of  the  approach,  the  vertical  tail  can  be  represented  as  a 
flexible  cantilevered  beam  with  rigidly  fixed  root  and  rigid  cross-section.  The  tail 
bending  Y  and  torsional  0  deflections  occur  about  an  elastic  axis  that  is  not  necessary 
coincident  with  the  inertial  axis.  If  spanwise  motion  is  ignored,  then  components  of  the 
velocity  of  die  beam  centre  of  gravity  can  be  described  in  the  inertial  frame  as 
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where  z  axis  of  die  frame  coincides  with  the  tail  elastic  axis  and  measures  the  vertical 
distance  from  the  fixed  support  along  die  tail  length  while  x  axis  is  placed  in  the 
vertical  tail  plane  and  xe  (z)  =  xCG  (z)  -  xe  (z)  is  the  distance  between  the  elastic  axis 
and  inertial  axis. 

Thus,  die  kinetic  mid  potential  energy  of  the  system  can  be  expressed  as 
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accordingly,  while  the  work  performed  by  the  aerodynamic  loading  can  be  expressed 
as 

L 

WA  =  J  {N(z,  t)Y  (z,  t)  +  M(z,  1)0  (z,  t))dz  (6) 

0 

where:  EI( z)  and  GJ{ z)  are  the  bending  and  torsional  stiffness  of  the  tail  section, 
g  is  the  acceleration  due  to  gravity, 

(3  is  the  angle  between  vertical  tail  and  delta  wing  plane, 
m( z)  is  the  mass  per  unit  length, 

/cg(z)  is  the  mass  moment  of  inertia  per  unit  length  about  the  centre  of  gravity, 
N( z)  is  the  normal  force  per  unit  length, 

M( z)  is  the  torsional  moment  per  unit  length. 

Introducing  Equations  (4)-(6)  into  Equation  (1),  the  following  differential  equations  of 
motion  can  be  obtained: 
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The  initial  assumption  of  small  motion  about  an  equilibrium  point  allows  for  tire  above 
equations  of  motion  to  be  presented  in  linear  form  as: 
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w/(z)xe  (z)g  cos  P  =  M(z,t) 

(10) 


For  a  flexible  vertical  tail  represented  as  a  cantilevered  beam  with  rigidly  fixed  root,  tire 
geometrical  and  natural  boundary  conditions  on  Y  and  0  are  given  as: 


17 


DSTO-RR-0256 


and 


mo 


8Y(  Q,Q  d2Y(L,i ) 
dt  ~  dz 2 


8_ 

8z 


EI(L) 


d2Y(L,t) 

dz2 


=  0 , 


0(0,0  = 


8Q(L,t) 

dz 


(11) 


(12) 


The  solution  of  Equations  (9)  and  (10)  with  the  above  boundary  conditions  provides 
the  total  response  of  a  dynamic  system  to  an  unsteady  aerodynamic  loading. 

Following  Strganac  [52],  the  system's  total  response  can  be  represented  in  terms  of  the 
static  and  dynamic  solutions  as 


and 


Y(z,t)  =  Ys(z)  +  Yd(z,t), 


e(z,o=0s(z)+omo- 


(13) 

(14) 


In  case  of  decaying  motions  the  solution  of  tire  Equations  (9)  and  (10)  will  converge  to 
tire  static  solution  and  may  not  require  unnecessary  computation  of  the  dynamic 
contributions.  Also,  for  the  case  of  large  sta+tic  deformations,  the  solutions  of  the 
dynamical  equations  about  the  static  deformations  will  conform  to  tire  original 
assumption  of  small  motions  about  an  equilibrium  points.  Therefore,  the  Equations  (9) 
and  (10)  can  be  presented  in  terms  of  the  static  and  dynamic  contributions  as 


m{z) 


d2Yd(z,t) 


82 


8z 2 


dt2 


Em 


+  w(z)*o  (z) 


d2ed(z,t)  d2 


dt2 


+ 


dz2 


Em 


dz2 


8z 2 


-  m(z)g cos  P  =  N s (z)  +  Nd(z,t), 


and 


, ,  mmo ,  ,  ^2Yd(z,t)  8  f „r,  ,50,(z)' 

4(z) - - - +  m(  z)x„(z) - — - —  GJ(z)  - — 

dt  dt"  dz\  dz  j 

8  (  <90  ( 7  /  ) N 

GJ(z) — -/w(z)gjfe(z)cosp  =  Ms(z)  +  Md(z,t). 


dz 


dz 


(15) 


(16) 


Here,  the  equations  for  tire  static  deformations  can  be  written  as 


d2 


8z~ 


Em 


dz2 


-m(z)gcosp  =  Ns(z), 


(17) 


o 

dz 


50,( z) 


GJ(z) 

\  dz 


-  z)x0  (z )g  cos  P  =MS  (z), 


(18) 
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and  the  dynamic  equations  for  small  motion  about  the  static  deformations  can  be 
presented  as 


d%M  ,  ,d2Qd(z,t)  ,  d2 


m{  z)- - - - +  w(z)x0  (z)  ,  . 

dt2  dt2  dz2 


EI(z) 


V 


8%(z,Q 

dz2 


h(z) 


ae,(M)  _ 527,(z,0  df_S0rf(z,^ 


w(z)xe(z)- 


<9/“  dz 


GJ(z) 

v  Sz  j 


=  Nd(z,t),  (19) 


=  Md(z,t).  (20) 


A  solution  of  the  above  equations  can  be  obtained  by  an  expansion  of  the  dependent 
variables  in  terms  of  the  natural  free  vibration  modes  of  the  system,  see 
Meirovitch  [63].  Here,  the  variables  Yd  and  can  be  represented  as 

r,(z,  0  =  i>,  (*)■?,  (0>  (21) 

/=1 

and 

Mz>0  =  Jjv  j(*)q  j(t),  (22) 

j=i+\ 


where  \p,  and  \|/j  are  modeshapes  satisfying  the  free-vibration  modes  of  bending  and 
torsion, 

q,  and  q2  are  generalised  coordinates  for  bending  and  torsion  modes, 

/  is  the  number  of  bending  modes, 

(J-I)  is  die  number  of  torsion  modes  selected  to  represent  the  solution. 

Substituting  these  expansions  into  Equations  (15)  and  (16)  and  using  the  Galerkin 
method  along  with  integration  by  parts  and  imposing  die  boundary  conditions. 
Equations  (11)  and  (12),  one  can  get  the  following  equations  for  die  generalised 
coordinates  q,  and  q,  in  matrix  form: 


1 

n 

l _ 

'4,  +  F/ 

1 

Is) 

Is) 

Is) 

1 _ 

s 

~  Mn 

Mn 

(  n2  X 
°  <1, 

dt 2 

_  7 

+ 

X 

1 - 

ri 

V 

M2l 

M  22 

d-qj 

l  dt2  ) 

d 

X 

k22\ 

fjj 

f 


V 


(23) 


(24) 


where 
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[K]is  the  stiffness  matrix, 

[M]  is  the  mass  matrix, 

(As )  is  tire  static  component  of  aerodynamic  loading  vector, 

(A(I  )  is  the  dynamic  component  of  aerodynamic  loading  vector, 
(F)  is  tire  static  loading  vector  due  to  gravity, 

and 


Lr  <3\|/  .  d~\\l  . 

Ku  =  f  £/(»)-“■  aV<fe 

for  /  =  1,2,...,/  and  j  =  1,2,...,/ 

k12  =  k2]=o 

Kn  = 

•  oz  oz 

for  /  =  1  + 1,...,  J  and  j  =  I + 

L 

M  I,  =  jdz 

0 

for  7  =  1,2,...,/  and  j  =  1,2,...,/ 

l 

Mn  =  M21  =  Jw(z)x0(z>|/,\|/Jr/z 

0 

for  /  =  1,2,...,/  and  j  =  /  + 1,...,  J 

L 

M22  =  J/0  (*>l'  A'  jdz 

0 

for  /  =  7  +  1,...,  J  and  j  =  /  +  1,...,J 

Ai  =  jy\>,Nx(z)dz 

0 

for  /  =  1,2,...,/ 

L 

A,j  =  jMJz)dz 

0 

for  j  =  I  + 

L 

A,  =  \\\fiNd(z,t)dz 

0 

for  /  =  1,2,...,/ 

/. 

A  =  jMd(z,t)dz 

0 

for  y'  =  /  +  l,...,J 

L 

F,  =  ^\i,m(z)gdz 

0 

for  ?'  =  1,2,...,/ 

i. 

Ft  =  |v|/  tm(z)gx0  (z)dz 

for  j  =  /  + 1,...,  J 

o 


which  are  similar  to  those  developed  by  Strganac  [52].  Note,  that  tire  above  equations  of 
motion  are  coupled  both  iner daily  and  aerodynamically,  as  the  sectional  center  of 
gravity  and  elastic  axis  are  not  geirerally  coincident  and  the  aerodynamic  loads  are 
motion  dependent. 
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4.  Unsteady  Aerodynamic  Model 

The  equations  of  motion  (23,  24)  describing  dynamic  response  of  the  system  require  an 
aerodynamic  model  that  is  able  to  predict  unsteady  buffet  loads  on  a  vertical  tail 
simultaneously  with  the  motion.  Although  Euler  or  Navier-Stokes  codes  can  be  applied 
for  the  simulation  of  unsteady  flows  over  delta  wings,  the  limitation  of  the  available 
computational  resources  calls  for  more  economical  methods  when  it  comes  to  such 
computationally  intensive  problems  as  time-accurate  prediction  of  empennage  buffet. 
Therefore,  a  vortex  approach  is  used  for  the  development  of  an  unsteady  aerodynamic 
model  that  provides  economy  in  computation  by  focusing  the  computational  effort  on 
the  areas  of  concentrated  vorticity. 


4.1  Overview  of  Vortex  Method 

Vortex  methods  simulate  the  unsteady  fluid  flows  under  the  assumption  of  nonlinear 
dynamics  of  vorticity,  so  fluid  flows  at  high  Reynolds  numbers  can  be  simulated  by 
regions  of  concentrated  vorticity  embedded  in  ir rotational  fluid.  The  evolution  of  vortex 
elements  can  be  tracked  numerically  in  a  Lagrangian  or  hybrid  Euler-Lagrangian 
reference  frame.  As  stated  by  the  theorems  of  Helmholtz  and  Kelvin,  the  inviscid  motion 
of  vorticity  in  these  regions  is  provided  by  the  local  fluid  velocity,  which  in  turn  is 
determined  kinematically  from  the  vorticity  field.  Thus,  it  is  very  convenient  to  consider 
inviscid  fluid  dynamics  in  terms  of  concentrated  vorticity  regions,  which  induce  motion 
on  each  other  as  an  alternative  to  pressure-velocity  considerations. 

Over  the  past  decades,  vortex  methods  have  been  developed  and  applied  as  a  reliable 
predictive  tool  for  a  wide  range  of  engineering  problems.  Their  area  of  application  has 
extended  fr  om  simple  inviscid  potential  flow  calculations  to  the  modelling  of  unsteady 
viscous  flows  by  simulating  the  three-dimensional  Navier-Stokes  equations.  This  ability 
of  vortex  methods  to  account  for  the  complex  features  of  fluid  dynamics  from  within  a 
single  computational  framework  makes  them  a  very  attractive  tool  for  a  wide  range  of 
engineering  problems. 

It  was  shown  by  Levinski  [64,  65],  that  the  vortex  model  is  able  to  predict  tire  onset  of 
leading  edge  vortex  breakdown  past  a  delta  wing,  twin- tail  configuration  and 
qualitatively  describe  its  unsteady  behaviour.  Also,  spatial  and  temporal  characteristics 
of  unsteady  buffet  pressures  arising  on  the  rigid  tail  were  found  to  be  in  qualitative 
agreement  with  available  experimental  data.  A  detailed  description  of  the  unsteady 
vortex  method  can  be  found  in  Levinski  [65]  and  only  a  brief  review  of  the  technique  is 
provided  here. 
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4.2  Governing  Equations 

Fluid  flows  are  governed  by  the  Navier-Stokes  equations,  which  for  incompressible 
viscous  fluid  flow  may  be  expressed  as 


—  +  (u-V)u  =  -^-+vV2u,  (25) 

dt  p 

V  •  u  =  0 ,  (26) 


where  u  is  a  flow  velocity, 
p  is  a  pressure, 
p  is  tire  fluid  density, 
v  is  a  flow  viscosity. 

In  homogeneous  fluids,  vorticity  to  is  defined  as 

to  =  V  x  u.  (27) 

Then  in  terms  of  vorticity.  Equation  (25)  can  be  written  as 

^  +  (u-V)co  =  (o-Vu  +  vV2co,  (28) 

which  is  the  vorticity  transport  equation  and  represents  simultaneous  convection  and 
diffusion  of  vorticity  in  the  flow  field  as  well  as  concentration  of  vorticity  due  to  vortex 
filament  stretching.  The  solution  can  be  obtained  using  the  widely  adopted  fractional  step 
method,  which  simulates  tire  convection  and  the  diffusion  processes  sequentially  rather 
than  simultaneously,  see  Leonard  [66],  This  can  be  achieved  by  solving  Euler's  equations 
by  the  vortex  method  and  using  the  proper  numerical  technique  to  model  tire  diffusion 
equation.  The  required  local  velocities  can  be  computed  as  the  solution  of  Poisson's 
equation  for  tire  velocity  field  in  terms  of  Biot-Savart  integration. 


4.2.1  Vorticity  Production  and  Convection 

In  a  real  flow,  the  well-known  viscous  boundary  layer,  caused  by  tire  presence  of 
viscosity,  develops  adjacent  to  tire  body  surface.  The  essence  of  the  majority  of  vortex 
methods  is  to  replace  this  boundary  layer  by  a  region  of  concentrated  vorticity,  which 
causes  a  reduction  in  the  fluid  velocity  from  its  local  value  to  zero  right  on  tire  body 
surface.  For  computational  simplicity  the  real  fluid  flow  can  be  represented  as  a  bulk  of 
irrotational  inviscid  outer  flow  separated  from  the  actual  body  surface  by  tire  thin 
bouirdary  shear  layer,  where  tire  effects  of  viscosity  and  turbulence  are  simulated. 

For  veiy  high  Reynolds  number  flows  tire  fluid  can  be  considered  as  inviscid  and  tire 
boundary  layer  cannot  be  developed  on  tire  body  surface.  In  such  flows  tire  fluid 
velocity  changes  discontinuously  from  zero  at  the  point  that  belongs  to  the  body 
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surface  up  to  the  value  of  local  flow  velocity  just  above  it.  In  this  case,  the  body  surface 
itself  must  be  represented  as  the  distributed  vorticity  sheet  a  of  strength  y  and  die 
boundary  condition  of  zero  normal  velocity  is  used  to  ensure  'no-leak'  velocity  through 
the  body  surface. 

It  can  be  shown  that  the  rate  of  vorticity  production  in  the  shear  layer  is  directly  related  to 
the  pressure  gradient  [66],  so  vorticity  creation  is  largely  attributable  to  the  dynamic 
behaviour  of  the  outer  flow.  Here,  a  statement  of  the  Neumann  boundary  condition  of 
'no-leak'  velocity  at  a  point  k  on  a  body  surface  a  may  be  expressed  by  the  following 
integral  equation 


-J-|g(K,M)r(n)rfcj„xn,  +/„(!!,, V„)  =  0,  (29) 

a 

where  nJC  is  a  unit  normal  vector  at  point  K  . 

p  is  a  variable  point  on  a  body  surface  c  , 

Q  is  a  velocity  singular  kernel, 

fn  is  a  function  of  surface  geometry  and  outer  flow  velocity. 

The  integral  equation  (29)  describes  the  process  of  vorticity  creation  on  the  solid 
boundaries  and  can  be  used  to  define  the  distribution  of  vorticity  over  the  body  surface.  It 
states  that  the  local  vorticity  density  y  (p)  must  be  such  that,  in  the  presence  of  an 

incident  flow  velocity  Vco,  the  resultant  flow  remains  parallel  to  the  body  surface. 

A  natural  consequence  of  tire  process  of  vorticity  production  and  convection  along  the 
body  surface  is  its  shedding  from  certain  separation  points  into  the  flow  field.  Although 
in  a  real  flow  the  viscous  boundary  layer  may  separate  spontaneously  under  various 
conditions,  the  assumption  of  a  fully  attached  boundary  layer  is  commonly  used  in 
inviscid  potential  flow  calculations.  In  such  an  approach  the  flow  is  allowed  to  separate 
only  at  a  finite  number  of  prescribed  separation  points,  such  as  tire  body  trailing  edge,  in 
order  to  satisfy  the  unsteady  Kutta-Joukowski  condition,  see  Belotserkovskii  and 
Nisht  [6  7\. 

The  equations  of  motion  of  separated  vortices  are  given  by  the  value  of  Are  velocity  field 
at  their  present  locations 


~  =  u(r,,0-  (30) 

at 

The  spatial  evolution  of  the  flowfield  can  be  obtained  by  integrating  the  above  equation 
using  the  Adams-Bashford  scheme,  which  is  recommended  for  vortex  methods  [66], 

The  required  local  velocities  can  be  computed  as  tire  solution  of  Poisson's  equation  for  tire 
velocity  field 
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V2u  =  -Vxoj.  (31) 

For  a  flow  field  without  interior  boundaries  and  with  the  fluid  at  rest  at  infinity,  the 
solution  to  Poisson's  equation  (31)  may  be  written  as  the  Biot-Savart  integral  equation, 
which  gives  a  velocity  distribution  induced  by  an  area  of  concentrated  vorticity,  see 
Lamb  [68],  Hence,  one  can  get  a  nonlinear  system  of  ordinary  differential  equations 
giving  tlie  temporal  evolution  of  the  vorticity  field. 


4.2.2  Discretisation  of  Vorticity  Regions 

According  to  the  vortex  method,  the  compact  vorticity  regions  can  be  confined  to  a 
finite  number  of  isolated  tubes  of  vorticity  or  vortex  filaments.  Similar  to  tire  point 
vortex  representation  in  two-dimensional  vortex  flows,  these  flows  in 
three-dimensional  coordinates  can  be  modelled  by  a  system  of  space  curves,  each  with 
zero  cross-sectional  area  and  constant  circulation. 

Although  tire  vortex  filaments  with  an  infinitesimal  core  is  tire  simplest  method  for 
vorticity  discretisation,  many  investigators  have  used  vortices  with  cores  of  constant  and 
finite  shape  or  vortex  'blobs',  see  for  example  Chorin  [69]  or  Kuwahara  &  Takami  [70].  In 
tire  vortex  blob  method,  the  vortex  core  has  a  constant  axisynrmetrical  shape  which  is 
defined  by  a  'cut-off  radius  and  a  'cut-off'  function  to  yield  a  description  of  vorticity 
distribution  within  tire  core.  Thus,  if  the  spatial  configuration  of  tire  vortex  filament  with  a 
blob-like  core  is  defined  by  space  curve  r,  (s,t) ,  then  its  vorticity  field  is  given  by 


.dr, 


w,(r,/)  =  r,  J/G [r - r, (s, /)] ds. 


(32) 


where  s  is  a  parameter  along  tire  curve, 
fa  is  a  cut-off  function, 
r,  is  a  vortex  filament  circulation  . 

Following  Leonard  [71],  the  velocity  field  induced  by  this  filament  in  air  unbounded 
domain,  with  no  interior  boundaries,  can  be  determined  as 


dr, 


u'(r>')=-iTrJ 


[r-r,0M)]x— ^(|r-r,.(s,0|/5,) 


ds 

r-r ,(s,t) 


-ds. 


(33) 


where  smoothing  function  c  is  determined  by  tire  choice  of  'cut-off  function,  hr  the 
current  study,  tire  curvilinear  vortex  filament  of  arbitrary  shape  is  represented  as  a 
combination  of  rectilinear  vortex  filaments  of  constant  strength. 
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The  use  of  distributed  vortex  cores  or  vortex  blobs  yields  more  realistic  vorticity 
distributions  and  bounded  induced  velocities  compared  to  the  vortex  filaments  with 
infinitesimal  cores  [69].  It  can  be  shown  that  the  vortex  blob  calculation  scheme  satisfies 
the  principles  of  conservation  of  circulation,  of  impulse  and  momentum  of  impulse  for  the 
system  of  vortices  with  the  same  core  sizes,  see  Leonard  [71]. 

Despite  certain  advantages  of  the  high-order  explicit  cut-off  functions,  a  simpler 
function  of  vorticity  distribution  within  the  core  as  proposed  by  Rosenhead  [72]  is 
adopted  in  the  present  study.  The  core  radius  is  assumed  to  be  time  invariant  and 
similar  for  all  the  vortex  elements.  Such  a  scheme  was  found  to  be  adequate  for 
accurate  representation  of  the  velocity  field  modelled  by  well-spaced  vortex  filaments 
with  low  rate  of  adjacent  cores  overlapping. 


4.3  Computational  Algorithm 

To  derive  a  solution  for  the  unsteady  aerodynamic  problem  the  general  approach 
needs  to  be  formulated  with  appropriate  initial  and  boundary  conditions  as  well  as 
conditions  of  spatial  and  temporal  discretisation  of  the  body  and  its  wake.  A  numerical 
solution  is  developed  in  the  following  sections  by  the  author  where  the  vortex  filament 
method  is  used  to  describe  inviscid,  incompressible  flow  aerodynamics.  The  solution  is 
derived  using  the  fractional  step  method,  which  provides  a  sequential  simulation  of  the 
physical  processes  of  vorticity  development  and  evolution  in  the  flow  field. 


4.3.1  Formulation  of  the  Unsteady  Problem 

Consider  a  lifting  body  a  ,  undergoing  an  arbitrary  motion  in  an  infinite  space  of 
homogeneous  incompressible  fluid.  As  the  governing  equations  of  motion  are  developed 
in  terms  of  a  translating,  but  non-rotating,  reference  frame,  then  it  is  appropriate  to  define 
the  unsteady  flow  problem  in  the  body-fixed  coordinate  system  that  is  aligned  with  the 
vertical  tail  root  and  has  linear  velocity  U0(f)  with  respect  to  tire  fixed  inertial  frame. 

In  the  general  case  the  body  motion  can  be  described  by  its  linear  velocity  U(r,  /)  and 
angular  velocity  Q(r,  t)  with  respect  to  the  body-fixed  reference  frame.  Therefore,  the 
resulting  velocity  at  a  point  on  a  body  surface  a  can  be  defined  as 

W(rc  ,/)  =  V(r,/)  +  U0(/)  +  U(rc ,/)  +  fi(r0  ,/)  x  ra ,  (34) 


where  V(r,t)  is  the  velocity  field,  induced  by  the  body  and  its  wake. 

Here,  the  unknown  velocity  field  V(r,f)  in  the  fluid  domain  must  satisfy  the 
Navier-Stokes  equation 
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+  (V(r ,/)  •  V)V(r,/)  =  -- Vp(r,t)  +vV2  V(r  ,t),  (35) 

of  p 


as  well  as  the  continuity  equation 


V  •  V(iy)  =  0.  (36) 

It  is  known  that  for  high  Reynolds  number  flows  convective  processes  in  the  wake 
dominate,  and  tire  influence  of  viscous  diffusion  is  minimal.  For  such  flows  with 
negligibly  small  viscosity,  tire  assumption  of  fully  inviscid  flow  can  be  justified  and 
adopted.  In  this  case,  tire  original  equation  (35)  in  the  inviscid  fluid  domain  can  be 
reduced  to  Euler's  equation  for  the  unknown  velocity  field  V(r,/) 


dV(r  t )  1 

-AjLi  +  (V(r,/) .  V)V(r,/)  =  —Vp(r,t).  (37) 

dt  p 


Tire  induced  velocity  field  V(r,/)can  be  assumed  to  be  irrotational  in  the  whole  fluid 
domain  except  tire  body  surface  a  and  its  trailing  wake  a  i.e., 

(  \ 


V  x  V(r,/)  =  0, 


r  g  o  u 


(38) 


V  J 


Then  for  tire  iirduced  velocity  field  a  velocity  potential  O  exists  such  drat 

V(r,/)  =  V<D(r,0.  (39) 

Therefore,  accordiirg  to  Equations  (36)  and  (39)  tire  unknown  induced  velocity  field  can  be 
found  as  a  solution  to  Laplace's  equation 


AO(r,0  =  0, 


r  £  a  u 


(40) 


Tire  fundamental  objective  of  such  a  problem  is  to  find  tire  velocity  potential  0(r,  t)  that 
satisfies  Laplace's  equation.  Under  tire  above  initial  and  boundary  conditions,  the 
problem  has  a  unique  solution,  which  gives  us  tire  unknown  velocity  potential  0(r,f)  in 
tire  fluid  domain  [68],  It  is  also  known  from  Lamb  [68],  that  tire  solution  for  Laplace's 
equation  can  be  obtained  by  a  suitable  choice  of  basic  singularities  such  as  doublets, 
sources  or  vortices.  In  fact,  no  one  among  these  singularities  has  a  special  advantage  over 
tire  others,  as  tire  potential  flow  can  be  modelled  correctly  by  air  appropriate  distribution 
of  either  sources  or  vortices.  This  is  guaranteed  by  tire  uniqueness  theorems  for  Laplace's 
equation.  But  tire  numerical  solution  can  be  formulated  at  different  levels  of  accuracy  and 
computational  efficiency.  Therefore,  it  is  particularly  advantageous  to  consider  the  flow 
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simulation  method,  based  on  vortex  singularities,  because  such  important  physical 
processes  as  flow  separation  from  the  body  surface  as  well  as  unsteady  wake 
development  and  their  interaction,  can  be  modelled  using  a  single  approach. 


4.3.2  Representation  of  the  Body  and  its  Wake 

Following  vortex  theory,  a  body  surface  and  its  wake  can  be  replaced  by  infinitely  thin 
layers  of  vorticity  or  vortex  sheets  consisting  of  spanwise  and  chordwise  vortex  lines  of 
variable  strength.  Here,  the  presence  of  spanwise  vorticity  in  a  body  vortex  system  is 
caused  by  a  continuous  change  of  tire  boundary  condition  at  the  wall  brought  about  by 
the  flow  unsteadiness.  The  unsteadiness  of  the  boundary  conditions  is  also  responsible  for 
the  presence  of  spanwise  vorticity  in  a  continuous  vortex  wake  separating  from  the  body 
surface.  In  the  process  of  spatial  discretisation  of  the  continuous  vorticity  layer,  the 
curvilinear  vortex  lines  of  variable  intensity  can  be  represented  as  a  system  of  rectilinear 
chordwise  and  spanwise  vortex  filaments  of  constant  strengths. 

In  the  current  numerical  scheme  the  boundary  layer  is  allowed  to  separate  only  at  a  finite 
number  of  separation  points,  assuming  potential  flow  modelling  for  the  rest  of  the  body. 
Here,  the  vorticity  diffusion  and  convection  activities  within  the  body  boundary  layer  are 
completely  ignored,  thereby  substantially  reducing  the  computational  effort.  It  is  assumed 
that  the  places  of  separation  are  known  a  priori  and  fixed  at  sharp  edges  of  the  body 
surface.  Such  an  approximation,  as  well  as  the  assumption  of  fully  attached  flow  on  tire 
rest  of  the  body  surface,  is  commonly  used  in  inviscid  flow  calculations.  Following  this 
simplified  technique  a  reasonably  accurate  prediction  of  overall  aerodynamic 
characteristic  may  be  obtained  especially  for  sharp  edged  bluff  bodies. 


4.3.3  Initial  and  Boundary  Conditions 

As  the  process  of  flow  development  will  be  considered  from  the  time  t  =  0,  it  is  assumed 
that  during  tire  period  of  time  t  <  0  the  fluid  in  the  domain  and  the  lifting  body  itself  are  in 
a  state  of  rest.  Then  for  t  >  0  both  the  flow  and  the  body  are  suddenly  set  in  motion. 
According  to  vortex  theory,  the  body  unsteady  motion  will  be  accompanied  by  tire 
development  of  its  trailing  wake  which  can  be  represented  as  a  combination  of  vortex 
wake  sheets  a  .  Therefore,  tire  unknown  velocity  field  V(r ,t)  in  Equation  (35)  must 
satisfy  tire  initial  condition 

V(r,0  =  0,  W(r,/)  =  0,  t  <  0,  (41) 

as  well  as  the  following  boundary  conditions: 

1.  For  iir viscid  flow  a  statement  of  the  Neumann  boundary  condition  of  'no-leak' 
velocity  at  any  point  on  a  body  surface  a  can  be  written  as 

W(ro ,/)  -n(r0  ,/)  =  0,  (42) 
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where  n(ra  ,/)  is  a  unit  normal  vector  to  die  body  surface  at  ra  . 

2.  At  an  infinite  distance  from  the  body  and  its  wake,  the  condition  of  induced 
velocities  decay  must  be  satisfied  as 

limV(r,O  =  0.  (43) 

r-^oc 

This  condition  is  satisfied  identically  through  the  Bio-Savart  Law. 

3.  Tire  Kutta-Joukowski  condition  of  finite  velocities  must  be  satisfied  at  the  lines  p  of 
trailing  wake  separation  such  tiiat 

W(rp,t)  *  oo.  (44) 


As  a  consequence  of  die  Kutta-Joukowski  condition,  a  vortex  wake  will  be  created  at 
die  points  of  dow  separation.  The  kinematic  condition  of  normal  velocity  continuity 


w, (r„  ,0 •  »,(r  ,0  =  W,(r  ,0-n/r  t), 


(45) 


as  well  as  pressure  continuity 


Pi(ro,’0  =  Pj(rOf,0,  (46) 

must  be  imposed  on  the  upper  i  and  lower  j  surfaces  of  the  vortex  wake  sheet  a  . 

Following  die  Kelvin-Helmholtz  dieorem,  diis  kinematic  condition  is  satisfied  by 
allowing  die  trailing  wake  to  move  with  the  local  fluid  velocity  while  preserving  its 
circulation,  see  Lamb  [68]. 


4.3.4  Calculation  of  die  Pressure  Distribution 


Having  obtained  die  induced  velocity  potential  0(r,i)  from  the  solution  of  Equation  (40), 
die  unsteady  velocity  field  W(r,/)  can  be  determined.  The  unsteady  pressure  field  in  die 
inertial  frame  can  be  calculated  using  die  Cauchy-Lagrange  equation  [67]  for  inviscid 
incompressible  dow  as 


Ptr.n=tjn 


ot  2 


(47) 


where  die  pressure  function  fp  (t)  can  be  determined  from  die  boundary  condition  at 
infinity.  Thus,  if  die  induced  velocities  decay  at  infinite  distance  from  the  body  and  its 
wake,  so  tiiat 


0(r,/)  — »  0,  |V0(r,/)|  — »  0, 


(48) 
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and  assuming  constant  pressure  at  infinity,  one  can  get 

Ar,0  =  P.~  (49) 

ot  2 

Therefore,  an  unsteady  pressure  coefficient  on  the  body  surface  a  can  be  defined  as 

f  w  (T»0  2  ,  (50) 

PUq/2  ot 

where  U0  -  is  the  flow  field  characteristic  velocity, 

w  =  W/  U0  -isa  dimensionless  flow  velocity, 

<f)  =  <J>  /  U06  -isa  dimensionless  velocity  potential. 

The  unsteady  pressure  distribution  can  be  found  using  the  potential  difference  A<|)  across 
the  vorticity  sheet  a  ,  which  is  defined  as 

(Xa )  =  <|>f.  (ra )  -  (r„ )  =  ,  (51) 

where  T0  is  the  circulation  of  a  vortex  element  representing  the  vortex  sheet  at  ro  ,  and 
indices  i  and;  refer  to  the  upper  and  lower  sides  of  the  vorticity  sheet  a  ,  respectively. 


Thus,  the  Equation  (50)  can  be  presented  as 
Cp(ro,0  =  l-w2(ra,0-2 


5[<j)(ro,0-<t)(r,,0]  -3<l>(r„/) 


where  rs  is  a  coordinate  vector  of  the  leading  edge  stagnation  point. 

Then,  according  to  the  Equation  (51) 

<l>,  (ra  (r,  ,0  =  W>7  (rc , t)  +  Acf) (rc , /)]  -  [<t>7  (r4 , /)  +  A(J> (r, , /)].  (53) 

Assuming  that  A<()(ro  )  =  ro  and  A<j)(r5 )  =  Fv ,  the  above  equation  can  be  differentiated 
as 

(r,,0]  =  ^[<t>7(ro,0-(t>7(rj,0]  +  -^(ro  -T,).  (54) 

ot  ot  ot 


According  to  the  boundary  conditions,  there  is  no  motion  allowed  inside  the  body 
surface,  so 

<t>.(rc,O-<Mrs,O  =  0,  (55) 
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giving 

■^[♦,(r„,0-+,(r„0]  =  -|(ro  -r,). 

at  at 


(56) 


Due  to  temporal  discretisation  of  tire  problem,  the  above  equation  can  be  written  as 


By  definition 


at 


C-n+r;-1). 

At 


a 

<Mr*>0 =<!>«  +  JWA> 


(57) 

(58) 


where  ([y  is  the  potential  of  free-stream  velocity  field  and 


3<t>,(r ,»0 


Bt 


J- 


dt 


-dx. 


(59) 


Haviirg  defined  all  tire  rurkirown  terms  from  Equation  (50),  a  distribution  of  unsteady 
pressure  coefficients  can  be  integrated  over  tire  whole  body  surface  to  obtain 
instantaneous  values  of  aerodynamic  forces  N(z,t)  and  moments  M(z,t )  required  for 
the  solution  of  the  equations  of  motion  (23),  (24). 


5.  Numerical  Solution  of  Aeroelastic  Equations 


Tire  aeroelastic  equations,  describing  dynamic  response  of  a  buffeting  tail,  and  an 
unsteady  aerodynamic  model  that  determines  motion-dependent  buffet  loads  have  been 
developed  in  non-dimensional  form  in  previous  sections.  After  introduction  of  the 
variable  ry  =  dqi{t)ldl ,  tire  governing  equations  of  motion  (23),  (24)  result  iir  a  set  of 
first-order  differential  equations,  which  can  be  presented  in  state  space  form  as 


o  -M44IJ  Ml  {4,1 

w  o  Jlwui  o  I’ 


(60) 
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see  Strganac  [52].  Here,  the  number  of  equations  is  twice  the  number  of  modes  chosen  to 
represent  the  dynamical  system. 

Solution  of  the  above  equations  in  the  time  domain  presents  a  certain  challenge  since 
the  aerodynamic  loads  depend  on  the  tail  motion,  yet  the  tail  motion  cannot  be 
determined  unless  the  aerodynamic  loads  are  known.  Following  Strganac,  an  iterative 
numerical  integration  scheme  that  accounts  for  the  interaction  between  the 
aerodynamic  loads  and  tail  dynamic  response  is  developed  which  determines  the 
motion  of  the  structure  and  the  motion-dependent  aerodynamic  loads  simultaneously. 
The  method  is  based  on  a  fourth-order  Hamming's  predictor-corrector  method,  see 
Carnahan,  et  at  [73],  An  advantage  of  the  predictor-corrector  method  is  that  it  does  not 
subdivide  a  time  step  of  integration,  as  do  other  numerical  integration  techniques,  such  as 
the  Runge-Kutta  method,  providing  substantial  economy  of  computations. 

Hamming's  scheme  requires  the  values  of  the  dependent  variables  at  the  current  and 
three  previous  time  steps  as  well  as  their  derivatives  at  the  current  and  two  previous  time 
steps.  Using  initial  conditions,  the  starting  values  of  the  dependent  variables  and  their 
derivatives  can  be  found  using  a  fourth  order  Runge-Kutta  method.  Having  this  data,  the 
values  of  state  variables  f:  for  the  next  time  step  t  + 1  can  be  predicted  by  a  fourth  order 
Milne  predictor  as 


fi'+l  =//"3  +yA/(2//'  -fr'+if;'-2) 


(61) 


At  tire  end  of  each  time  step,  die  local  truncation  error  e'  is  estimated  and  used  to 
modify  the  predicted  value  of  the  state  variables  as 


~  117 

f'+ 1  -  fl+l  +llZ.p' 

J  7,0  J  7  ,0  T  g  ci  ’ 


(62) 


assuming  that  the  value  of  the  local  truncation  error  does  not  change  significantly  on 
successive  intervals  [73].  Then  Hamming's  corrector  equation  is  solved  iteratively 
using  the  successive  substitution  algorithm  at  each  iteration  step  j  to  obtain  the 
updated  state  variables 


f!j' -ik/r -/rJ+3A'fe-i+2/"  ~/r\  j =u . * 


8 


(63) 


until  convergence  of  all  the  state  variables  is  achieved  after  k  iterations. 

Once  convergence  of  the  state  variables  and  loads  has  been  achieved,  the  local 
truncation  error  for  the  corrector  equation  is  estimated  as 
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and  die  final  values  of  die  state  variables  are  determined  similarly  to  Milne's  method 
as 


fr  =  fi 


/+i 

,k 


/+ 1 


(65) 


For  each  time  step,  integration  of  die  equation  of  motion  by  Hamming's  mediod  begins 
widi  die  computation  of  die  tiuid  dow  by  solving  Equation  (37)  and  convection  of  wake 
vorticity  to  die  new  force-free  position  from  the  position  generated  at  the  end  of  the 
last  time  step.  A  dowchart  detailing  the  general  approach  for  computation  of  the 
dynamic  solution  is  presented  in  Figure  1. 

Here,  the  new  state  variables  are  predicted  using  aerodynamic  loads  computed  at  the 
last  time  step  and  applied  to  die  structural  dynamics  model.  This  new  wake  geometry 
and  new  state  variables  are  used  to  compute  the  aerodynamic  loads.  Then,  the  state 
variables  are  corrected,  new  tail  geometry  is  computed,  the  aerodynamic  loads  are 
updated,  and  convergence  of  both  die  state  variables  and  aerodynamic  loads  is 
checked.  If  convergence  is  not  achieved,  the  state  variables  are  corrected  again  and  new 
aerodynamic  loads  are  determined.  This  process  is  repeated  until  convergence  of  the 
loads  and  tail  motion  is  achieved.  Here,  the  new  state  variables  are  predicted  for  the 
first  integration  pass  and  are  corrected  for  all  subsequent  passes  while  maintaining  the 
dowfield  in  a  fixed  position. 

Next,  the  wake  vorticity  is  re-convected  from  its  position  generated  at  the  end  of  the 
last  time  step  using  dow  conditions  based  on  the  initial  and  corrected  values  of  die 
state  variables.  This  updated  position  of  die  wake  vorticity  is  checked  for  convergence 
witii  its  previous  position  and  is  convected  by  an  iterative  process  to  a  new  converged 
position  while  keeping  the  state  variables  fixed.  The  iteration  continues  with  die  use  of 
die  new  converged  dowfield  conditions  for  updating  die  state  variables  until  their 
convergence  with  die  aerodynamic  loads  is  achieved  during  the  fixed  dowfield 
computations  and  then  again  proceeds  witii  the  iterative  convection  of  wake.  This 
two-stage  iterative  scheme  is  aimed  at  providing  full  convergence  of  die  wake  position, 
state  variables  and  die  aerodynamic  loads  at  the  end  of  each  time  step  prior  to 
advancing  to  die  next  time  step. 

Despite  tiiis  integration  scheme  for  die  aeroelastic  equations  being  able  to  provide 
ratiier  accurate  convergence  of  motion  and  loads,  its  excessive  computational 
requirements  were  found  to  be  prohibitive  for  buffet  computations  when  long 
computed  time  histories  are  required  to  resolve  time-dependent  data.  Thus,  a 
simplified  algorithm  was  investigated  for  integration  of  the  equations  of  motion,  where 
the  geometry  of  the  wake  remains  the  same  as  generated  at  the  first  iteration  pass  and 
is  not  updated  each  time  die  state  variables  are  determined,  see  Figure  2.  When  die 
new  tail  geometry  is  computed,  the  aerodynamic  loads  are  updated,  and  convergence 
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of  both  the  state  variables  and  aerodynamic  loads  is  checked.  If  convergence  is  not 
achieved,  the  state  variables  are  corrected  and  new  aerodynamic  loads  are  determined. 
Here,  only  unsteady  aerodynamic  loads  are  updated  each  time  the  state  variables  are 
predicted  or  corrected.  This  process  is  repeated  until  convergence  of  the  loads  and  tail 
motion  is  achieved  while  maintaining  the  flowfield  in  a  fixed  position. 

It  was  found  that  the  results  of  the  lower-order  integration  scheme  for  some  test  cases 
show  no  noticeable  difference  over  the  more  accurate  higher-order  scheme  while 
providing  substantial  economy  of  computation.  Therefore,  this  simplified  integration 
scheme  was  adopted  and  subsequently  used  for  modelling  the  vertical  tail  buffet  in  this 
study. 


5.1  Temporal  and  Spatial  Discretisation 

Note,  that  the  vortex  algorithm  for  simulation  of  unsteady  separated  flow  employs  an 
iterative  time  marching  scheme  in  which  one  system  of  vortex  elements  is  shed  from  the 
body  surface  at  every  time  step.  The  accuracy  of  such  a  model  depends  on  a  correct 
balance  between  the  scales  of  spatial  and  temporal  discretisation  of  the  unsteady  problem. 
For  high  Reynolds  number  flows  with  dominating  influence  of  convection  over  diffusion, 
an  appropriate  spatial  discretisation  of  the  body  must  correspond  to  the  scale  of  vorticity 
convective  motion.  According  to  vorticity  inviscid  dynamics,  an  average  convective 
displacement  A/c  of  die  vorticity  in  a  shear  layer  can  be  defined  as 

(66) 


where  U0  is  a  free-stream  velocity  and  A/c  is  time  step  based  on  vorticity  convective 
motion.  A  link  between  the  convection  displacement  of  vorticity  A/c,  governed  by  inviscid 
vorticity  dynamics,  and  the  corresponding  scale  As  of  the  body  surface  discretisation  may 
be  defined  as 


k 


(67) 


where  A:  is  in  the  range  of  [0...1].  Here,  the  scale  of  the  body  surface  discretisation 
depends  on  the  number  of  vortex  elements  nv  =  lc/As  approximating  the  body  surface 
where  lc  is  the  body  characteristic  length  equal  to  the  length  of  the  delta  wing  root 
chord. 


Here,  the  scale  of  temporal  discretisation  must  be  chosen  such  as  to  satisfy  numerical 
stability  of  the  integration  scheme  for  the  structural  dynamics  part  of  the  problem  in 
order  to  resolve  the  highest  order  modes  selected  to  represent  the  dynamic  system. 
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6.  Computation  of  Vertical  Tail  Dynamic  Response 

While  an  ultimate  aim  of  the  work  is  to  investigate  the  F/A-18  empennage  buffet,  it 
was  shown  by  Levinski  [65]  that  the  vertical  tail  buffet  problem  could  be  simulated 
efficiently  using  a  generic  delta  wing,  twin  vertical-tail  configuration.  This  simple 
configuration  contains  all  tire  pertinent  physics  involved  in  the  development  and  burst 
of  a  leading-edge  vortex  and  its  subsequent  interaction  with  the  vertical  tails.  Here,  the 
delta  wing  is  used  to  generate  the  leading-edge  vortices,  which  tend  to  break  down 
over  the  wing  resulting  in  highly  turbulent  swirling  flow.  This  energetic, 
vortex-breakdown  flow  then  impinges  on  the  vertical  tails  producing  unsteady, 
unbalanced  buffet  loads  on  the  tail  surfaces  and  causing  their  dynamic  response 
(buffeting). 

The  dyiramic  aeroelastic  computational  model  described  in  the  previous  sections  is 
employed  to  simulate  the  development  of  the  unsteady  buffet  loads  produced  by  the 
vortex  breakdown  flow  as  well  as  to  predict  the  vertical  tail  dynamic  response.  In  this 
way,  tire  major  characteristics  of  the  vertical  tail  buffet  can  be  investigated  without 
complications  associated  with  simulation  of  the  flow  over  the  complete  F/A-18  aircraft. 


6.1  Computational  Model 

Tire  same  delta  wing,  twin  vertical-tail  configuration  as  used  in  the  previous  buffet 
study  of  Levinski  [65]  is  employed  to  computationally  simulate  the  vertical  tail 
dynamic  respoirse.  It  consists  of  a  sharp-edged,  76-degree  leading  edge  sweep  delta 
wing  and  swept-back  F/A-18  twin  tails,  see  Figure  3  and  Figure  4.  Both  tire  delta  wing 
and  twin  tails  are  of  zero  thickness.  Each  of  tire  tails  is  of  aspect  ratio  1.2  with  the  root 
chord  length  of  0.4  and  a  tip  chord  length  of  0.159,  based  on  the  delta  wing  root  chord. 
The  tails  are  cantilevered  on  tire  upper  surface  of  a  trailing  edge  extension  of  the  delta 
wing  and  have  a  sweepback  angle  of  35  degrees  for  a  quarter-chord  spanwise  line.  The 
tip  of  each  tail  is  inclined  outboard  such  that  dihedral  angle  between  the  two  tails  is  40 
degrees.  The  spanwise  separation  distance  between  tire  tails  at  tire  root  is  50%  of  the 
delta  wing  span. 

The  numerical  simulation  procedure  requires  a  discretisation  of  the  computational 
model,  thus  a  total  of  2986  closed  vortex  rings,  or  panels,  are  used  for  spatial 
discretisation  of  tire  delta  wing,  trailing  edge  extension  and  twin  vertical  tails.  A  unit 
aspect  ratio  delta  wing  is  modelled  using  2178  quadrilateral  vortex  panels.  Tire  vortex 
system  of  tire  trailing  edge  extension  continues  the  discretisation  pattern  of  tire  delta 
wing  and  is  represented  by  416  rectangular  panels  while  196  trapezoidal  panels  model 
each  of  the  vertical  tails. 

The  'no-penetration'  boundary  condition  is  imposed  on  all  of  the  lifting  surfaces.  The 
unsteady  wake  is  modelled  as  a  vortex  sheet,  composed  of  vortex  panels,  which 
emanates  from  tire  sharp  edges  of  tire  delta  wing,  trailing  edge  extension  and  vertical 
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tails.  Only  symmetrical  flow  cases  are  considered  in  the  present  buffet  study.  In  all  the 
test  cases,  the  condition  of  flow  symmetry  is  enforced  during  the  computational 
simulation  of  the  flow  over  the  delta  wing,  twin-tail  configuration. 

Each  of  the  elastic  vertical  tails  of  the  computational  model  is  treated  as  a  swept  back 
beam  which  is  allowed  to  oscillate  in  bending  and  torsion  modes  representing  the  first 
bending  and  torsion  modes  of  the  F /  A-18  vertical  tails.  In  the  present  study  of  coupled 
bending-torsion  response,  the  distance  between  the  elastic  axis  and  the  inertia  axis  for 
each  of  the  tails  is  set  equal  to  0.04  based  on  delta  wing  root  chord. 


6.2  Test  Conditions 

Previous  sub-scale  and  full-scale  experimental  investigations  of  the  F/A-18  vertical  tail 
buffet  as  well  as  flight  trials  have  shown  that  the  frequency  characteristics  and 
intensity  of  the  tail  buffet  vary  primarily  as  a  function  of  angle  of  attack  and  dynamic 
pressure.  Analysis  of  the  dynamic  flight  test  results  indicated  that  vertical  tail  first 
bending  (15.8Hz)  and  first  torsion  (45.7Hz)  modes  have  the  most  significant  impact  on 
the  fatigue  life  of  the  empennage  structure  so  these  two  modes  are  considered  as 
fatigue  critical  modes  of  vibration.  The  flight  conditions  at  which  peak  response  levels 
occur  for  the  most  damaging  vertical  tail  modes  are  found  to  be  at  20-39  degrees  angle 
of  attack  and  175-400  psf  dynamic  pressure  [74],  Based  on  the  critical  dynamic  response 
region,  the  test  conditions  for  the  dynamic  aeroelastic  model  were  selected  to  vary 
between  20  and  40  degrees  angle  of  attack  at  250  psf  dynamic  pressure  and  altitude  of 
10,000  ft  as  these  conditions  are  most  important  in  vertical  tail  buffet  studies. 

During  the  simulations,  the  fluid  was  started  impulsively  and  computation  was  carried 
out  for  300  time  steps  with  a  dimensionless  time  increment  of  At  -  0.012  that  allows  the 
vortex  wake  to  travel  8.6  tail  root  chord  lengths  at  the  rate  of  the  free-stream  velocity. 
At  the  selected  test  condition  of  10,000  ft  altitude  and  250  psf  dynamic  pressure,  the 
total  computed  time  interval  corresponds  to  0.143  seconds  of  real  flight  time.  This 
computational  time  span  is  deemed  to  be  adequate  to  obtain  a  fully  developed  vortex 
wake  allowing  the  unsteady  buffet  loads  to  be  considered  as  random  and  stationary  in 
a  statistical  sense. 

However,  at  the  above  test  conditions  the  dynamically  scaled  elastic  tail  would 
complete  only  two  oscillations  in  the  first  bending  mode,  and  that  is  not  sufficient  for  a 
reliable  statistical  analysis  of  its  dynamic  response  especially  when  taking  into  account 
the  initial  transition  period  of  vortex  wake  development.  Therefore,  stiffness  and 
inertial  properties  of  the  elastic  tail  are  selected  such  that  during  the  computational 
time  span  the  vertical  tail  is  expected  to  complete  more  than  four  oscillations  in  the  first 
bending  mode  and  about  eight  oscillations  in  the  first  torsion  mode. 
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6.3  Vortex  Wake  Structure 

Particle  traces  of  unsteady  wake  development  over  the  delta  wing,  vertical-tail 
configuration  at  different  angles  of  attack  after  300  time  steps  are  presented  in  Figure  5, 
Figure  7  and  Figure  9  for  the  rigid  tail  computations  obtained  in  [65]  and  in  Figure  11, 
Figure  13  and  Figure  15  for  the  flexible  tail  computations.  Here,  only  the  leading  edge 
vortex  wakes  are  visualised  for  clarity.  Comparison  of  the  results  shows  a  similarity  of 
the  wake  structures  for  all  the  test  cases.  As  can  be  seen  in  all  tire  test  cases,  the  initially 
stable  leading  edge  vortex  experiences  perturbations  and  breaks  down  prior  to 
impinging  on  tire  vertical  tails  and  results  in  a  highly  disturbed  wake,  which  convects 
downstream  undergoing  gradual  expansion. 

The  characteristic  feature  of  the  wake  structure  for  all  the  test  cases  is  the  presence  of 
periodic  vortex  clusters  in  tire  severely  degenerated  vortex  wake.  The  wake  periodicity 
is  observed  in  the  results  obtained  during  both  the  rigid  and  flexible  tail  computations. 
The  presence  of  coherent  structures  in  the  wake  indicates  an  existence  of  a  dominant 
frequency  in  the  process  of  wake-structure  interaction.  It  appears  that  dynamic 
response  of  the  flexible  tails  affects  the  process  of  formation  of  the  vortex  clusters,  as 
their  periodicity  is  different  compared  with  the  rigid  tail  results. 

Details  of  the  leading  edge  flow  development  are  provided  in  Figure  6,  Figure  8  and 
Figure  10  for  rigid  tail  computation  and  in  Figure  12,  Figure  14  and  Figure  16  for 
flexible  tail  computation.  Here,  the  leading  edge  vortex  is  visualised  by  particle  traces, 
which  are  plotted  for  several  consecutive  time  steps.  As  one  can  see,  for  all  the  test 
cases  the  vortex  model  reproduces  the  leading  edge  vortex  roll-up,  which  is  typical  of 
the  flow  over  delta  wings.  This  leading  edge  vortex  expands  in  a  continuous  manner 
while  convecting  downstream.  The  vortex  core  is  initially  intact  and  stable  as  indicated 
by  die  coinciding  particle  traces.  However,  at  some  point  downstream,  the  growing 
scatter  of  tire  particle  traces  reveals  the  onset  of  instability  and  perturbation  of  die 
vortex  cores,  followed  by  sudden  enlargement  of  the  vortex  core  size. 

For  both  die  rigid  and  dexible  tail  calculations,  die  increase  of  die  vortex  core  size  is 
followed  by  development  of  highly  diffuse  and  disorganised  vortical  tiow,  which 
indicates  an  appearance  of  die  vortex  breakdown.  However,  in  die  case  of  die  flexible 
vertical  tail,  die  regions  of  vortex  wake  appear  to  be  larger  and  more  diffuse  as  a  result 
of  interaction  of  wake  vorticity  widi  die  deflecting  vertical  tails. 

It  can  be  seen  for  all  the  test  cases  that  highly  diffuse  vortex  flow  convects  downstream 
and  partially  covers  die  region  inboard  of  die  twin  vertical  tails.  The  twin  tails  cut 
dirough  die  leading  edge  vortices,  forming  two  vortical  tiows  inside  the  region 
between  the  twin  tails  and  smaller  vortical  flows  outside  each  tail.  As  in  die  case  of 
rigid  tails,  die  vortical  tiows  inside  the  region  are  larger  but  weaker  than  tiiose  vortices 
drat  develop  outside  the  tips  of  die  twin  tails.  It  is  observed  from  die  plots  of  vortex 
wake  structure  diat  interaction  of  die  leading  edge  vortex  flow  with  the  flexible  vertical 
tail  causes  more  degeneration  of  the  vortex  wake  compared  with  the  rigid  tail  case,  as 
indicated  by  highly  scattered  particle  traces. 
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The  severity  of  leading  edge  vortex  burst  becomes  more  pronounced  at  the  higher 
angle  of  attack.  For  example,  at  40  degrees  angle  of  attack,  the  computed  leading  edge 
vortex  is  larger  in  size  and  its  deterioration  into  an  unsteady  wake-like  flow  starts 
earlier  compared  to  the  lower  angle  of  attack  cases,  see  Figure  16.  As  indicated  by  a 
progressing  scatter  in  particle  traces,  this  perturbation  in  the  well-defined  leading  edge 
vortex  structure  produces  more  severe  degeneration  of  the  vortex  wake. 

Insignificant  scatter  of  the  particle  traces  at  the  20-degree  angle  of  attack  in  Figure  12 
shows  only  slight  perturbation  of  the  leading  edge  vortex  when  it  approaches  the 
trailing  edge  region.  It  indicates  that  in  this  case  the  wake  has  not  deteriorated  to  the 
degree  found  at  the  higher  angles  of  attack  and  actually  remains  stable  all  the  way 
along  the  wing.  This  result  is  consistent  with  the  result  of  rigid  tail  computation 
presented  in  Figure  6  and  suggests  that  for  this  angle  of  attack  tire  location  of  vortex 
breakdown  is  very  close  to  the  trailing  edge. 

Analysis  of  results  indicates  that  for  the  angle  of  attack  of  20  degrees  and  higher,  the 
predicted  burst  location  is  forward  of  the  wing  trailing  edge  and  moves  upstream  with 
the  increase  in  die  angle  of  attack.  In  all  the  test  cases,  the  breakdown  flow  is  located 
inside  the  region  between  die  twin  tails.  It  can  be  concluded  that  the  increase  in  the 
angle  of  attack  changes  the  shape  of  the  vortex  breakdown  dow  and  the  area  occupied 
by  the  breakdown  dow  becomes  larger  in  size. 

Additional  investigation  into  the  structure  of  the  leading  edge  vortex  breakdown  is 
performed  using  a  three-dimensional  view  of  the  iso-surfaces  of  time-averaged 
pressure  coefficient  at  30  degrees  angle  of  attack,  see  Figure  17  and  Figure  18,  for  rigid 
and  dexible  tails,  respectively.  As  indicated  by  the  iso-surfaces  in  Figure  18,  the  centre 
of  the  vortex  breakdown  dow  is  located  inboard  of  the  vertical  tail  which  is  consistent 
with  the  rigid  tail  results  and  predictions  of  others,  see  [59,  65].  Similar  to  the  rigid  tail 
case,  the  area  of  vortex  breakdown  grows  in  size  as  it  approaches  the  vertical  tail  and 
die  decelerated  dow  starts  to  recover  just  after  it  passes  the  wing's  hailing  edge. 
However,  it  is  evident  that  the  dedections  of  the  vertical  tail  change  the  location  and 
shape  of  die  vortex  breakdown  dow.  The  area  of  dow  deceleration  is  larger  in  size  and 
placed  slightly  upstream  indicating  that  the  leading  edge  vortex  is  experiencing  more 
intense  breakdown  compared  witii  die  rigid  tail  case.  Also,  the  duid-structure 
interaction  becomes  more  intense  near  the  tail  tips  where  dedections  of  die  tail  are  die 
largest,  and  that  results  in  further  deceleration  of  vortex  breakdown  dow  in  die  upper 
tail  region. 


6.4  Vertical  Tail  Buffet  Loading 

The  magnitude  and  frequency  content  of  instantaneous  pressure  ductuations  on  die 
tail  surface  are  the  most  important  characteristics  of  tail  buffet.  Analysis  of 
time-dependent  random  buffet  data  was  performed  using  statistical  data  reduction 
techniques  in  the  time  and  frequency  domains,  see,  for  example,  Pettit,  et  al.  [75], 
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The  differential  (buffet)  pressure  was  estimated  from  time  histories  of  surface  pressure 
fluctuations  as  the  difference  between  the  inner  surface  pressure  values  and  outer 
surface  pressure  values  on  tire  tail.  The  values  of  unsteady  differential  pressures  were 
integrated  over  tire  vertical  tail  surface  to  give  components  of  aerodynamic  load 
coefficients.  The  unsteady  pressure  and  loads  coefficients  were  reduced  to 
root-mean-square  (RMS)  and  power  spectral  density  (PSD)  forms.  The  peak  power  and 
the  dominant  frequency  of  the  unsteady  buffet  loads  were  determined  from  their  PSD 
plots. 

The  root-mean-square  (RMS)  value  of  zero-mean  differential  pressure  AP'(r)  at  the 
vertical  tail  provides  a  measure  of  magnitude  of  the  time-averaged  fluctuations  of  the 
unsteady  net  pressure  across  the  tail.  It  was  determined  by  averaging  the  differential 
pressure  histoiy  over  a  specified  elapsed  time  te  -  /s  as  follows 

AP'!(r)  =  — !—  [(/>,, (r.O -  C(r,<))!  dl  =— !—  )&P'(r,r)dt,  (68) 

where  /*„(r,/)and  Pom(r,t)  are  zero-mean  inner  and  outer  pressure  values  on  the  tail. 

The  averaging  is  performed  from  the  moment  the  wake  is  considered  as  fully  developed 
until  the  end  of  computation  so  that  the  dimensionless  times  for  start  ts  and  end  te  of 
sampling  are  selected  as  1.4  and  8.6  units,  respectively. 

The  RMS  differential  pressure  coefficient  is  then  defined  as 

^’(r)  -  AP'(r)/ qx,  (69) 

where  qx  is  the  free-stream  dynamic  pressure. 

The  frequency  content  of  instantaneous  pressure  fluctuations  was  estimated  using  the 
power  spectral  density  function  of  tire  unsteady  component  of  tire  differential  pressure 
coefficient.  Tire  differential  pressure  time  histories  from  each  test  condition  were 
converted  into  the  frequency  domain  using  Fast  Fourier  Transform  (FFT)  techniques. 
The  pressure  time  histories  were  divided  into  blocks  and  tire  Hanning  window  was 
applied  to  reduce  bandwidth  leakage.  The  average  PSD  functions  were  obtained  by 
averaging  tire  fast  Fourier  transforms  of  each  of  tire  blocks  with  50%  overlap  in  order  to 
increase  statistical  confidence.  The  dominant  frequencies  of  tire  buffet  pressure  were 
identified  from  the  power  spectral  density  plots  with  dimensionless  frequency 
resolution  of  0.73  based  on  the  tail's  root  chord. 
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6.4.1  Magnitude  of  Buffet  Pressures 

It  is  common  in  vertical  tail  buffet  studies  to  examine  and  compare  characteristics  of 
surface  pressure  fluctuations  measured  at  the  45%  chord  and  60%  span  location  on  the 
vertical  tail.  It  provides  a  common  base  for  comparison  of  the  results  of  various 
experimental  and  numerical  tail  buffet  studies.  A  computed  time  history  of  the 
differential  pressure  fluctuations  at  45%  chord  and  60%  span  location  obtained  on  the 
flexible  vertical  tail  at  30  degrees  angle  of  attack  is  presented  in  Figure  19  and 
accompanied  by  computational  results  obtained  on  a  rigid  tail.  As  one  can  see,  for  both 
the  rigid  and  flexible  tails,  the  differential  pressure  coefficients  fluctuate  substantially 
and  exhibit  oscillatory  motion  with  large  random  fluctuations  about  a  local  mean 
value.  However,  the  magnitude  of  differential  pressure  fluctuations  on  the  flexible 
vertical  tail  is  lower  compared  with  the  rigid  tail  case.  This  reduction  of  magnitude  of 
buffet  pressure  fluctuations  measured  on  the  flexible  tail  is  consistent  with  the  trends 
obtained  from  experimental  and  flight  test  data  in  Figure  20,  [30].  Here,  the  full-scale 
wind  tunnel  test  of  F/A-18  tail  buffet  [75]  was  performed  at  low  dynamic  pressure  so 
the  tail  does  not  deflect  and  it  could  be  considered  as  almost  rigid  in  contrast  to  flight 
condition  where  the  tail's  dynamic  response  is  significant. 

As  one  can  see  in  Figure  20,  the  flight  and  full-scale  data  follows  the  same  trend  as 
angle  of  attack  increases,  reaching  a  maximum  value  of  buffet  pressure  near  30  degrees 
angle  of  attack.  However,  the  flight  data  are  lower  in  magnitude  than  the  full-scale  test 
data  through  all  the  angle  of  attack  range,  with  the  largest  difference  of  up  to  25%  near 
30  degrees  angle  of  attack,  which  is  explained  by  'elastic  relief  effects'. 

Variation  of  RMS  differential  pressure  fluctuations  obtained  at  45%  chord  and  60% 
span  location  on  rigid  and  flexible  tails,  and  a  full-scale  wind  tunnel  test  of  F/A-18  tail 
buffet  [75]  at  different  angles  of  attack  is  presented  in  Figure  21.  As  one  can  see,  the 
RMS  buffet  pressure  fluctuations  obtained  on  the  flexible  tail  are  lower  than  the  rigid 
tail  values  for  all  the  angles  of  attack  but  higher  than  full-scale  test  data.  Despite  the 
computed  RMS  buffet  pressures  not  being  expected  to  closely  match  measured  data 
due  to  simplified  geometry  of  the  computational  model,  the  overall  reduction  of  buffet 
pressure  values  representing  the  effect  of  tail  flexibility  matches  experimentally 
observed  trends. 


6.4.2  Spectral  Content  of  Buffet  Pressures 

Spectral  content  of  the  differential  pressure  fluctuations  at  30  degrees  angle  of  attack  at 
the  45%  chord,  60%  span  location  is  presented  in  Figure  22  for  rigid  and  flexible  tail 
models.  One  can  see  that  on  the  rigid  tail  the  differential  pressure  contains  energy  over 
a  relatively  narrow  frequency  band  with  centre  frequency  of  2.9  corresponding  to  the 
dominant  frequency  of  vortex  breakdown  flow.  This  differential  pressure  power  peak 
occurs  at  about  the  same  frequency  values  across  tire  whole  surface  of  the  rigid  tail. 
However,  power  spectral  densities  of  considerably  lower  magnitudes  are  predicted  for 
the  flexible  tail  model  where  two  distinct  peaks  are  now  aligned  with  tire  frequencies  of 
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die  tail  bending  and  torsion  modes.  The  presence  of  pressure  fluctuations  near  the  tail 
natural  frequencies  is  the  result  of  coupled  flow-structure  interaction  where  elastic  tail 
response  modifies  the  differential  pressure  fluctuations  initially  caused  by  turbulent 
vortex  breakdown  alone.  These  changes  in  magnitude  and  dominant  frequencies  of 
power  spectral  density  of  differential  pressure  between  the  rigid  and  flexible  tails 
found  at  the  45%  chord,  60%  span  location  appear  to  be  representative  for  the  majority 
of  other  locations  on  the  tail. 

Random  fluctuations  of  the  buffet  pressures  are  caused  by  interaction  of  the  highly 
turbulent  leading  edge  vortex  flow  with  the  vortex  system  of  the  twin  vertical  tails.  The 
leading  edge  vortex  breakdown  causes  flow  deceleration  and  tire  leading  edge  vorticity 
tends  to  concentrate  in  the  region  near  the  twin  tails,  gradually  building  up  a  vortex 
cloud.  The  process  of  formation  and  destruction  of  the  vortex  clouds  interacting  with 
the  tail  vortex  system  drives  the  oscillatory  behaviour  of  the  buffet  loads.  These 
drifting  clouds  of  scattered  vorticity  are  also  subjected  to  random-like  motion  caused 
by  vortices  bouncing  back  off  the  tail  surface  that  also  imposes  a  certain  degree  of 
randomness  to  the  surface  pressure  distribution. 

The  periodicity  of  the  vortex  cloud  formation  determines  the  dominant  frequencies  of 
buffeting  flow.  However,  this  is  further  affected  by  the  dynamic  response  of  the 
flexible  vertical  tails.  The  vertical  tail  deflections  change  the  location  and  shape  of  the 
vortex  breakdown  flow  that  modifies  pressure  distribution  on  the  tails,  as  well  as  its 
frequency  content,  by  the  appearance  of  motion-driven  pressure  fluctuations.  The 
resulting  buffet  pressure  magnitudes  and  frequency  content  depend  on  intensity  and 
coupling  of  both  the  unsteady  aerodynamic  loads  arising  from  tine  leading  edge  vortex 
burst  and  those  caused  by  tine  vibration  of  tine  tail. 

Note,  that  the  vortex  model  overpredicts  the  characteristic  frequency  of  the  pressure 
fluctuations  on  the  rigid  tail  as  tine  dominant  frequency  of  about  0.5  was  identified  in 
most  of  the  sub-scale  and  full-scale  experiments  on  F/A-18  tail  buffet ,  see  [12,  75].  This 
characteristic  frequency  is  a  property  of  the  flowfield  of  the  burst  LEX  vortex  and 
further  development  of  tine  vortex  model  is  required  to  rectify  this  discrepancy. 


6.4.3  RMS  Pressure  Distribution 

Contour  plots  of  tine  calculated  RMS  differential  pressure  fluctuations  on  a  flexible 
vertical  tail  are  presented  in  Figure  23  through  Figure  27  for  five  values  of  angle  of 
attack  between  20  and  40  degrees.  They  are  accompanied  by  rigid  tail  results  to 
investigate  variation  of  buffet  pressure  distribution  at  different  test  conditions. 

For  both  the  rigid  and  flexible  models,  tine  magnitude  of  RMS  differential  pressure  in 
most  of  tine  locations  increases  with  the  increase  of  angle  of  attack  due  to  an  increase  in 
the  strength  of  tine  leading  edge  vortex  before  its  subsequent  burst  as  well  as  the 
upstream  movement  of  tine  location  of  vortex  breakdown. 
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On  the  flexible  tail,  the  pressure  fluctuations  are  higher  at  the  leading  edge  as  this  area 
is  the  closest  to  the  breakdown  location,  while  the  pressure  fluctuations  tend  to 
decrease  towards  the  trailing  edge.  However,  on  a  rigid  tail,  the  decrease  in  differential 
pressure  towards  the  trailing  edge  is  rather  uniform  while  on  a  flexible  tail  the  pattern 
is  more  complex  and  irregular.  On  the  flexible  tail,  high  RMS  differential  pressure 
values  are  found  not  only  near  the  leading  edge  close  to  the  tail's  root,  but  also  at  the 
tail  tip  and  near  the  upper  part  of  the  trailing  edge.  The  largest  gradients  of  the  RMS 
differential  pressure  values  are  still  detected  in  the  leading  edge  area.  Although  the 
general  trend  of  decrease  of  buffet  pressure  RMS  towards  the  trailing  edge  is  present 
for  most  of  the  test  cases,  the  pressure  field  is  obviously  affected  by  the  tail  dynamic 
response,  especially  at  those  parts  of  the  tail  surface  where  deflections  are  the  largest. 

In  general,  the  patterns  of  the  intensity  of  the  unsteady  pressure  fluctuations  on  rigid 
and  flexible  tails  have  more  similarity  for  angles  of  attack  of  30  degrees  and  higher, 
with  some  considerable  differences  appearing  for  the  lower  values  of  incidence.  For 
example,  at  20  and  25  degrees  angle  of  attack,  the  magnitude  of  the  RMS  differential 
pressure  on  the  rigid  tail  gradually  decreases  towards  the  trailing  edge  in  the 
chordwise  direction  and  towards  the  tail  tip  in  spanwise  direction.  Here,  the  position 
of  the  centre  of  the  vortex  breakdown  flow  accounts  for  the  pattern  of  the  pressure 
fluctuations  where  the  lower  half  of  the  tails  experience  larger  RMS  pressure  levels 
than  those  on  the  rest  of  the  tail.  However,  dynamic  response  of  the  flexible  tail  affects 
tire  flowfield  around  the  tail  and  this  considerably  changes  the  pattern  of  buffet 
pressure  distribution  resulting  in  lower  RMS  pressure  values  than  those  obtained  on 
the  rigid  tail  model. 

It  should  be  noted  that  an  improved  vorticity  convection  algorithm  was  incorporated 
into  the  flexible  tail  model  following  computation  of  the  rigid  tail  pressures.  It  was 
found  that  higher  RMS  differential  pressure  values  in  this  area  on  the  rigid  tail  model 
were  caused  by  the  development  of  a  local  stagnation  zone  and  overproduction  of 
vorticity  at  the  junction  of  vertical  tail,  and  delta  wing  and  its  trailing  edge  extension 
vortex  systems.  This  shortcoming  was  rectified  in  the  vorticity  convection  algorithm  in 
the  flexible  tail  model  and  this  partially  explains  the  difference  between  the  rigid  and 
flexible  tail  results  in  the  area  near  the  leading  edge  near  the  tail  root. 

The  results  indicate  that  an  increase  in  the  angle  of  attack  intensified  tail  buffet  by 
increasing  the  RMS  levels  of  differential  pressure  fluctuations  and  expanded  the  region 
of  peak  loading  from  the  root  of  the  tail  toward  the  leading  edge  and  centre  of  the  tail. 
Thus,  at  angles  of  attack  of  30  degrees  and  higher,  the  distribution  of  the  differential 
pressures  on  the  rigid  tail  is  more  even  in  the  spanwise  direction  and  decreases 
towards  the  trailing  edge  mainly  in  the  chordwise  direction.  The  results  for  the  flexible 
tail  show  similar  trends  with  the  major  difference  in  the  upper  trailing  edge  area  where 
deflections  are  large  but  there  is  little  interaction  with  vortex  breakdown  flow  as  it 
passes  over  the  upper  part  of  tire  tail. 
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6.4.4  Vertical  Tail  Buffet  Loads 

At  each  time  step,  the  zero-mean  differential  pressures  were  integrated  over  the  tail 
surface  to  obtain  the  instantaneous  values  of  normal  force  and  root  bending  moment 
coefficients.  The  time  history  of  the  zero-mean  unsteady  normal  force  coefficient  was 
determined  as 


CW)  =  — "t- t)-P„(r„0M=-l-XCaF(r„0Al, 

r  A  /■'  i  Af  , 


(70) 


where  Ap  is  the  total  planform  area  of  the  tail.  A,  is  the  area  of  the  z-th  panel  on  the  tail 
and  tire  instantaneous  pressure  differential  CA/)  is  evaluated  at  the  centre  of  the  z-th 
panel.  Here,  tire  surface  pressure  over  tire  panel  area  was  assumed  to  be  constant  and 
equal  to  tire  pressure  computed  at  the  centre  of  the  panel,  which  is  defined  as  the  point 
of  intersection  of  two  lines  connecting  midpoints  of  the  opposite  sides  of  the  parrel. 

The  time  histoiy  of  zero-mean  root  bending  moment  was  evaluated  in  a  similar 
manner,  assuming  that  the  resultant  force  over  the  i-tlr  panel  is  acting  at  the  panel's 
centre 


CM  (0  =  - 7 - X^» (*•,  >0  -  Pm, (r  > 0U,1,  =  — —  X CA P (X„t)A,l , ,  (71) 

q,  AFcF  ,  Afcf  , 


where  Cp  is  the  tail  root  chord  and  /,  is  the  distance  from  the  tail  root  to  the  centre  of  the 
z-th  panel. 

Finally,  tire  zero-mean  RMS  values  of  fluctuating  normal  force  and  root  bending 
momexrt  coefficients,  which  provide  a  measure  of  tire  average  fluctuations  of  the 
unsteady  buffet  loads  at  tire  vertical  tail  are  determined  as 


and 


(72) 


(73) 


respectively. 

Tire  variation  of  RMS  normal  force  coefficients  for  rigid  and  flexible  tails  as  a  function 
of  angle  of  attack  is  presented  in  Figure  28,  and  compared  with  full-scale  experimental 
data  from  [75],  Tire  variation  of  RMS  root  bending  moment  coefficients  at  different 
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angles  of  attack  is  also  presented  in  Figure  29  and  compared  with  the  results  of 
full-scale  tail  buffet  test  [75]. 

Comparison  of  the  computed  RMS  normal  forces  and  bending  moments  in  Figure  28 
and  Figure  29  shows  that  flexible  tail  model  produces  higher  values  of  both  the  normal 
force  and  bending  moment  coefficients  compared  with  rigid  tail  model  despite 
predicting  lower  RMS  differential  pressure  values  for  most  of  the  tail  surface.  This  is 
because  differential  buffet  pressures  generated  by  the  flexible  tail  model  have  higher 
values  of  cross-correlation  over  the  tail  surface  as  the  flowfield  near  the  tail  is  modified 
by  tail's  motion. 

In  general,  bodi  the  flexible  and  rigid  tail  models  produce  values  of  RMS  buffet  loads 
which  are  reasonably  close  to  those  measured  during  full-scale  test.  Note  that  despite 
higher  magnitudes  of  RMS  buffet  loads  generated  by  the  flexible  tail  model,  their 
variation  with  angle  of  attack  closely  matches  the  trends  observed  in  the  rigid  tail 
results.  Initially,  both  the  computed  and  measured  values  of  RMS  buffet  loads  increase 
with  the  angle  of  attack.  This  is  caused  by  the  increase  of  the  RMS  differential  pressures 
for  most  of  the  locations  on  the  tail,  due  to  the  increase  of  leading  edge  vortex  strength 
at  higher  angles  of  attack.  However,  both  the  rigid  and  flexible  models  fail  to  follow  the 
experimentally  obtained  trend  for  angles  of  attack  of  30  degrees  and  higher  where  the 
values  of  experimentally  measured  RMS  buffet  loads  start  to  decline  after  reaching  a 
peak  value  at  32  degrees  of  angle  of  attack.  Here,  predicted  RMS  buffet  loads  for  both 
rigid  and  flexible  tail  experience  a  local  minimum  at  30  degrees  and  continue  to 
increase  at  higher  angles  of  attack.  Further  numerical  and  experimental  investigation  is 
required  to  determine  if  such  a  behaviour  of  predicted  RMS  buffet  loads  is  either 
specific  to  the  computational  model  geometry  or  to  the  predictive  capabilities  of  the 
vortex  model  at  higher  angles  of  attack. 

Power  spectral  density  functions  of  unsteady  normal  force  and  root  bending  moment 
coefficients  are  determined  from  their  time  histories  in  a  manner  similar  to  the  buffet 
pressure  PSD.  Variation  in  the  spectral  content  of  unsteady  normal  force  coefficient 
with  angle  of  attack  is  presented  in  Figure  30  and  Figure  31  for  rigid  and  flexible  tails, 
respectively.  Variation  in  the  spectral  content  of  unsteady  root  bending  moment 
coefficient  with  angle  of  attack  is  also  presented  in  Figure  32  and  Figure  33  for  both  the 
rigid  and  flexible  tails. 

Investigation  of  the  PSD  plots  for  the  rigid  and  flexible  tails  indicates  that  the  peaks  of 
the  power  spectral  density  of  normal  force  coefficient  can  be  detected  for  all  die  angles 
of  attack,  and  dominant  frequencies  at  which  die  normal  force  is  exerted  on  die  vertical 
tail  tend  to  decrease  at  higher  angles  of  attack.  It  can  be  noted  that  there  is  a  similarity 
of  variation  of  dominant  frequencies  with  angle  of  attack  between  die  rigid  and  dexible 
tail  cases,  however  spectral  content  of  normal  force  coefficient  on  the  flexible  tad  seems 
to  be  affected  by  its  dynamic  response  in  bending  and  torsion  modes. 

Here,  for  the  flexible  tad  at  20  degrees  angle  of  attack,  a  main  peak  is  generated  at  the 
reduced  frequency  of  3.6  whde  a  noticeable  peak  is  also  present  near  0.73,  which  is  the 
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frequency  of  tire  tail  bending  mode.  At  30  degrees  angle  of  attack,  a  strong  single  peak 
of  load  power  is  aligned  with  die  frequency  of  1.9  diat  is  close  to  the  frequency  of  the 
tail  torsion  mode.  For  bodi  the  rigid  and  flexible  tails  at  40  degrees  angle  of  attack,  the 
stronger  of  the  two  peaks  is  generated  at  the  frequency  of  the  tail  bending  mode  and 
the  other  is  at  die  frequency  of  3.6.  A  number  of  minor  peaks  can  also  be  identified  in 
the  PSD  plots  but  dieir  frequencies  are  well  above  the  vertical  tail  modal  frequencies 
and  diey  are  not  considered  as  critical  for  the  present  study. 

Variation  in  die  spectral  content  of  unsteady  root  bending  moment  coefficient  with 
angle  of  attack  in  Figure  32  and  Figure  33  experiences  similar  trends,  however  at  40 
degrees  angle  of  attack,  a  much  stronger  peak  of  bending  moment  power  is  generated 
at  die  frequency  of  3.6  which  becomes  comparable  in  magnitude  to  the  other  major 
peak  that  is  aligned  with  the  frequency  of  the  tail  bending  mode. 


6.5  Vertical  Tail  Dynamic  Aeroelastic  Response 

Computational  time  histories  of  generalised  coordinates  q  for  bending  and  torsion 
modes  of  vibration  for  die  vertical  tail  are  presented  in  Figure  34  tiirough  Figure  38  for 
five  different  angles  of  attack  in  20  to  40  degrees  range.  The  time  histories  show  near 
harmonic  variations  of  generalised  coordinates  for  both  bending  and  torsion  responses 
for  most  of  die  test  cases.  However,  it  can  be  noticed  that  periodicity  of  dynamic 
response  was  not  fully  developed  due  to  insufficient  length  of  computed  time  histories. 
This  is  most  evident  for  die  bending  mode  due  to  its  lower  natural  frequency,  but  can 
also  be  seen  for  die  torsional  mode  having  intermittent  oscillations  of  higher  than 
average  amplitude.  A  certain  degree  of  irregularity  in  die  time  histories  of  generalised 
coordinates  for  bending  and  torsion  can  be  explained  by  aerodynamic  as  well  as 
structural  coupling  between  die  modes. 

It  is  known  from  test  data  that  for  high  angle  of  attack  manoeuvring  conditions  the 
modal  response  of  die  vertical  tails  is  vulnerable  to  changes  from  one  dominant  mode 
to  anodier,  see  Triplett  [76],  Measured  variation  of  relative  magnitudes  of  peak  PSD 
associated  with  the  first  two  modes  of  the  F/  A-18  tail  is  presented  in  Figure  39.  A  rapid 
decline  in  relative  magnitude  of  the  torsion  mode  after  25  degrees,  dominance  of  die 
bending  mode  by  35  degrees  and  their  crossover  at  about  30  degrees  angle  of  attack, 
are  die  most  interesting  features  of  the  modal  response  variation. 

In  the  present  study,  tire  lengtii  of  computed  time  histories  does  not  allow  for  accurate 
assessment  of  PSD  peaks  values.  Thus,  die  variation  of  participation  of  die  bending 
and  torsion  modes  in  die  total  dynamic  response  of  die  tail  at  different  angles  of  attack 
is  evaluated  by  calculating  generalised  coordinate  RMS  values  for  each  of  die  modes, 
see  Figure  40.  The  results  indicate  that,  in  general,  die  predicted  variation  in  modal 
response  follows  die  experimentally  observed  trends  witii  die  torsion  mode  peaking 
around  25  degrees  angle  of  attack  and  bending  mode  peaking  at  35  degrees  angle  of 
attack.  Despite  some  scatter  in  the  computed  data,  bodi  the  increase  in  relative 
magnitude  of  the  bending  mode  and  decline  of  die  torsion  mode  at  higher  angles  of 
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attack  are  reproduced  correctly.  Note,  that  the  predicted  decline  of  relative  magnitude 
of  bending  mode  after  35  degrees  angles  of  attack  is  also  in  agreement  with 
experimental  data. 

Computed  time  histories  of  tail  forward  and  aft  tip  deflections  are  shown  in  Figure  41 
through  Figure  45  for  five  values  of  angle  of  attack.  Here,  the  variation  of  relative 
magnitudes  of  the  bending  and  torsion  modes  determines  the  overall  tip  response 
where  die  contribution  of  the  torsion  mode  is  only  noticeable  at  25  degrees  of  attack  or 
higher  while  the  bending  mode  dominates  the  overall  tip  response  for  all  the  test 
conditions.  The  difference  between  the  amplitudes  of  tips  deflection  is  the  result  of 
coupling  between  the  bending  and  torsion  modes  and  the  presence  of  static  aeroelastic 
deformations  of  the  tail.  As  one  can  see,  the  tail  aft  tip  experiences  larger  total 
deflections  compared  with  the  forward  tip. 


7.  Conclusions 

The  development  of  a  computational  model  for  the  prediction  of  the  unsteady 
aeroelastic  behavior  of  a  flexible  tail  under  buffet-induced  loads  has  been  described. 
The  multidisciplinary  problem  of  tail  buffeting  is  solved  in  a  time  domain  using  an 
unsteady  vortex  model  for  the  prediction  of  aerodynamic  loads  and  coupled  aeroelastic 
equations  for  the  bending  and  torsional  deflections  of  the  tail  which  are  resolved  using 
the  Galerkin  method.  The  set  of  ordinary  differential  equations  governing  the  motion 
of  the  dynamic  system  is  integrated  by  a  predictor-corrector  algorithm. 

A  dynamic  aeroelastic  analysis  of  empennage  buffet  is  performed  for  a  generic  delta 
wing,  twin  vertical-tail  configuration  at  high  angles  of  attack.  The  computational 
model  consists  of  a  76-degree  sharp-edged  delta  wing  along  with  F/A-18  twin  vertical 
tails,  where  only  the  twin  vertical  tails  are  flexible  and  modelled  using  a  linear  modal 
method. 

The  aim  of  the  study  is  to  simulate  the  coupled  fluid-structure  interaction  and  present 
time-accurate  integration  of  the  empennage  dynamic  aeroelastic  response  to  die 
unsteady  pressure  loads  caused  by  the  vortex  breakdown  present  at  the  high  angle  of 
attack  conditions.  Results  include  time  histories  of  unsteady  buffet  loads  and  dynamic 
response  of  die  tail.  Comparisons  are  made  with  the  rigid  tail  differential  pressures  and 
loads,  as  well  as  available  test  data. 

It  was  shown  that  the  numerical  solution  is  able  to  predict  onset  of  the  vortex  burst  and 
qualitatively  describe  the  characteristics  of  unsteady  buffet  loads  on  the  vertical  tail.  It 
is  found  that  inclusion  of  die  tail  flexibility  affects  the  shape  and  location  of  vortex 
breakdown  and  results  in  degeneration  of  the  leading  edge  vortex  structure  into  a 
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larger  arid  more  diffuse  flow  region  where  the  location  of  vortex  breakdown  tends  to 
move  forward  compared  to  the  rigid  tail  case. 

It  is  found  that  tire  magnitude  of  differential  pressure  fluctuations  on  a  flexible  vertical 
tail  is  lower  compared  to  tire  rigid  tail  case  and  this  is  consistent  with  experimental  and 
flight  test  data.  Despite  both  the  rigid  and  flexible  tail  models  over-predicting  the 
measured  data,  the  overall  reduction  of  buffet  pressure  values  representing  the  effect 
of  tail  flexibility  matches  experimentally  observed  trends. 

A  comparison  of  computational  results  with  available  wind-tunnel  and  flight  test  data 
of  F/A-18  tail  buffet  revealed  that  the  dynamic  aeroelastic  model  performed 
reasonably  well  in  simulating  the  spatial  and  temporal  characteristics  of  the  buffet 
loads.  Computational  analysis  of  the  flexible  tail  showed  a  difference  in  buffet  loading 
when  compared  to  the  rigid  tail  model.  Here,  both  the  tail  bending  and  torsion  modes 
had  a  noticeable  contribution  to  tire  buffet  load  magnitude  and  frequency  content  and 
this  emphasises  the  need  to  account  for  tail  flexibility  during  computations  of  the 
empennage  buffet. 

Investigation  of  the  power  spectral  density  plots  for  the  buffet  loads  indicates  that  the 
dominant  frequencies  at  which  the  normal  force  and  root  bending  moment  are  exerted 
on  the  vertical  tail  tend  to  decrease  at  higher  angles  of  attack  and  this  compares  well 
with  available  experimental  results. 

It  is  also  noted  that  the  predicted  variation  in  modal  response  of  the  flexible  vertical  tail 
follows  the  experimentally  observed  trends.  Despite  some  scatter  in  the  computed 
data,  both  die  increase  in  relative  magnitude  of  the  bending  mode  and  decline  of  the 
torsion  mode  at  higher  angles  of  attack  are  reproduced  correctly. 

Although  the  aeroelastic  model  captured  many  of  die  experimentally  observed  trends, 
analysis  of  die  power  spectral  density  of  die  predicted  differential  pressures  on  the  tad 
indicated  diat  die  characteristic  frequency  of  vortex  breakdown  is  higher  than  the 
measured  values.  Predictive  capabdities  of  the  vortex  model  should  be  refined  in  order 
to  correctly  predict  variation  of  buffet  loads  and  dynamic  response  of  the  tail  at 
different  angles  of  attack.  This  is  left  for  future  studies. 

The  present  study  demonstrates  die  use  of  an  unsteady  vortex  model  coupled  with  a 
structural  dynamics  model  for  time-domain  aeroelastic  analysis  of  aircraft  empennage 
buffeting.  The  metiiod  is  capable  of  predicting  major  unsteady  features  of  vortex 
breakdown  induced  buffet  loads  that  are  required  for  characterizing  empennage  buffet 
and  die  resulting  structural  response.  The  present  research  provides  advances  in  die 
predictive  capabdity  and  our  understanding  of  aircraft  empennage  buffet. 
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Figure  1  General  integration  scheme  for  dynamic  solution 
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Figure  2  Simplified  integration  scheme  for  dynamic  solution 


Figure  3  General  view  of  delta  wing ,  twin  vertical-tail  configuration 


Figure  4  F/A-18  Vertical  Tail  Geometry 
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Figure  5  Top  view  of  vortex  wake  structure  at  20  degrees  angle  of  attack  for  rigid  tail 
computation 


Figure  6  Particle  traces  of  the  leading  edge  vortex  core  at  20  degrees  angle  of  attack  for  rigid 
tail  computation 
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Figure  7  Top  view  of  vortex  wake  structure  at  30  degrees  angle  of  attack  for  rigid  tail 
computation 
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Figure  8  Particle  traces  of  the  leading  edge  vortex  core  at  30  degrees  angle  of  attack  for  rigid 
tail  computation 
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Figure  9  Top  view  of  vortex  wake  structure  at  40  degrees  angle  of  attack  for  rigid  tail 
computation 


Figure  10  Particle  traces  of  the  leading  edge  vortex  core  at  40  degrees  angle  of  attack  for  rigid 
tail  computation 
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Figure  11  Top  vieiv  of  vortex  umke  structure  at  20  degrees  angle  of  attack  for  flexible  tail 
computation 


Figure  12  Particle  traces  of  the  leading  edge  vortex  core  at  20  degrees  angle  of  attack  for 
flexible  tail  computation 
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Figure  13  Top  view  of  vortex  ivake  structure  at  30  degrees  angle  of  attack  for  flexible  tail 
computation 


Figure  14  Particle  traces  of  the  leading  edge  vortex  core  at  30  degrees  angle  of  attack  for 
flexible  tail  computation 
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Figure  15  Top  view  of  vortex  wake  structure  at  40  degrees  angle  of  attack  for  flexible  tail 
computation 


Figure  16  Particle  traces  of  the  leading  edge  vortex  core  at  40  degrees  angle  of  attack  for 
flexible  tail  computation 


60 


DSTO-RR-0256 


Figure  17  Iso-surface  of  time-averaged  pressure  coefficient  at  30  degrees  angle  of  attack  for 
rigid  tail  computation 


Figure  18  Iso-surface  of  time-averaged  pressure  coefficient  at  30  degrees  angle  of  attack  for 
flexible  tail  computation 
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Figure  19  Time  history  of  differential  pressure  coefficient  at  45%  chord  and  60%  span  location 


Angle  of  Attack  (degrees) 


Figure  20  Comparison  of  RMS  differential  pressure  coefficient  values  at  45%  chord  and  60% 

span  location  on  the  tail  for  full-scale  wind  tunnel  test  and  in  flight  (Meyn  et  al.  [30]) 
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Figure  21  Comparison  of  RMS  differential  pressure  values  at  various  angles  of  attack  at  60% 
span,  45%  chord  location  on  the  tail 
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Figure  22  Power  spectral  density  of  differential  pressure  fluctuations  at  45%  chord  and  60% 
span  location  on  tlte  tail 
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Figure  28  Comparison  of  RMS  normal  force  coefficient  values  at  various  angles  of  attack 
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Figure  29  Comparison  of  RMS  bending  moment  coefficient  values  at  various  angles  of  attack 


Figure  31  Variation  of  flexible  tail  normal  force  power  spectral  density  with  angle  of  attack 
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Figure  32  Variation  of  rigid  tail  bending  moment  power  spectral  density  with  angle  of  attack 
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Figure  33  Variation  of  flexible  tail  bending  moment  power  spectml  density  with  angle  of  attack 
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Figure  34  History  of  generalised  coordinates  at  20  degrees  angle  of  attack 


Figure  35  History  of  generalised  coordinates  at  25  degrees  angle  of  attack 
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Figure  36  History  of  generalised  coordinates  at  30  degrees  angle  of  attack 
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Figure  37  History  of  generalised  coordinates  at  35  degrees  angle  of  attack 
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Figure  38  History  of  generalised  coordinates  at  40  degrees  angle  of  attack 
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Figure  39  Measured  variation  of  flexible  vertical  tail  dynamic  response  xvith  angle  of  attack 
(Triplett  [76]) 
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Figure  40  Variation  of  relative  magnitude  of  bending  and  torsion  modes  ivith  angle  of  attack 
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Figure  41  Tail  tips  deflection  history  at  20  degrees  angle  of  attack 
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Figure  42  Tail  tips  deflection  history  at  25  degrees  angle  of  attack 
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Figure  43  Tail  tips  deflection  history  at  30  degrees  angle  of  attack 
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Figure  44  Tail  tips  deflection  history  at  35  degrees  angle  of  attack 
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Figure  45  Tail  tips  deflection  history  at  40  degrees  angle  of  attack 
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